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Abstract 

We review properties of g-orthogonal polynomials, related to their orthogonality, du- 
ality and connection with the theory of symmetric (self-adjoint) operators, represented 
by a Jacobi matrix. In particular, we show how one can naturally interpret the duality 
of families of polynomials, orthogonal on countable sets of points. In order to obtain 
orthogonality relations for dual sets of polynomials, we propose to use two symmetric 
(self- adjoint) operators, representable (in some distinct bases) by Jacobi matrices. To 
illustrate applications of this approach, we apply it to several pairs of dual families of q- 
polynomials, orthogonal on countable sets, from the g-Askey scheme. For each such pair, 
the corresponding operators, representable by Jacobi matrices, are explicitly given. These 
operators are employed in order to find explicitly sets of points, on which the polynomials 
are orthogonal, and orthogonality relations for them. 
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"We mathematicians are particularly fond of duality theorems; translating 
mathematical statements from one category to another often gives us new 
and unexpected insight", M.Harris, "Postmodern at an Early Age", Notices 
of the American Mathematical Society, Vol.50, No. 7, p. 792, 2003. 

1. INTRODUCTION 

It is well known that each family Pnix), n = 0, 1, 2, • • •, of orthogonal polynomials of one 
variable corresponds to the determinate or indeterminate moment problem. If a polynomial 
family corresponds to the determinate moment problem, then there exists only one positive 
orthogonality measure n for these polynomials and they constitute a complete orthogonal set 
in the Hilbert space L^(^). If a family corresponds to the indeterminate moment problem, then 
there exists infinitely many orthogonality measures fi for these polynomials and these measures 
are divided into two parts: extremal measures and non-extremal measures. If a measure /i 
is extremal, then the corresponding set of polynomials constitute a complete orthogonal set 
in the Hilbert space L^(/i). If a measure /i is not extremal, then the corresponding family of 
polynomials is not complete in the Hilbert space i^(^) (see [ST]). 

It is also well known that there exists a close relation of the theory of orthogonal polyno- 
mials with the theory of symmetric (self-adjoint) operators, representable by a Jacobi matrix. 
The point is that with each family of orthogonal polynomials one can associate a closed sym- 
metric (or self-adjoint) operator A, representable by a Jacobi matrix. If the corresponding 
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moment problem is indeterminate, then the operator A is not self-adjoint and it has infinitely 
many self-adjoint extensions. If the operator A has a physical meaning, then these self-adjoint 

extensions are especially important. These extensions correspond to extremal orthogonality 
measures for the same set of polynomials and can be constructed by means of these measures 
(see, for example, [Ber], Chapter VII, and [Sim]). If the family of polynomials corresponds 
to the determinate moment problem, then the corresponding operator A is self-adjoint and 
its spectrum is determined by an orthogonality relation for the polynomials. Moreover, the 
spectral measure for the operator A is constructed by means of the orthogonality measure for 
the corresponding polynomials (see [Ber], Chapter VII). 

In section 2 we briefly review the relations between the theory of orthogonal polynomials, 
the theory of operators, representable by a Jacobi matrix, and the theory of moment problem. 
This information is basic for the exposition in further sections. In section 2 wc also discuss how 
one can naturally extend the conventional notion of duality to sets of polynomials, orthogonal 
on countable sets of points. 

In order to find orthogonality measures for dual sets of polynomials, we use two symmetric 
(or self-adjoint) operators, representable (with respect to different bases) by Jacobi matrices. 
This approach is applied to several sets of dual g-orthogonal polynomials from the Askey 
scheme. 

Pairs of operators (A, I), employed for studying some sets of g'-orthogonal polynomials 
and their duals, belong to the discrete series representations of the quantum algebra Uq{sui^i) 
(see, for example, [AKl] and [AK6]). However, in order to facilitate case of comprehending to 
a larger number of readers we have not exploited this deep algebraic fact; that is, we exhibit 
explicit forms of these operators without using the representation theory of the quantum 
algebra Uq{sni^i). These pairs of operators are, in fact, a generalization of Leonard pairs, 
introduced by P. Terwilliger [Terl] (for the definition and references see section 3). 

When one considers dual sets of g-polynomials, orthogonal on countable sets of points, 
then one member of these sets corresponds to the determinate moment problem and another 
to the indeterminate moment problem. One of the two operators (A, I) (that is, the opera- 
tor A) for a given dual pair of sets of g-orthogonal polynomials corresponds to a three-term 
recurrence relation for the set of polynomials, which corresponds to the determinate moment 
problem. This operator is bounded and self-adjoint; moreover, it has the discrete spectrum. 
We diagonalize this self-adjoint bounded operator and find its spectrum with the aid of the 
second operator /, which corresponds to a g-difference equation for the same set of polyno- 
mials. An explicit form of all eigenfunctions for the operator A is found for each dual set of 
polynomials, considered by us. They are expressed in terms of g-polynomials, which belong to 
the set, associated with the determinate moment problem. Since the spectrum of A is simple, 
its eigenfunctions form an orthogonal basis in the Hilbert space. One can normalize this ba- 
sis. This normalization is effected by means of the second operator / from the corresponding 
pair. As a result of this normalization, two orthonormal bases in the Hilbert space emerge: 
the canonical (or the initial) basis and the basis of eigenfunctions of the operator A. They 
are interrelated by a unitary matrix U, whose entries Umn are explicitly expressed in terms 
of polynomials Pm{x), corresponding to the determinate moment problem. Since the matrix 
U is unitary (and in fact it is real in our case), there are two orthogonality relations for its 
elements, namely 



The first relation expresses the orthogonality relation for the polynomials Pm{x), which cor- 
respond to the determinate moment problem. So, the orthogonality of U yields an algebraic 
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proof of orthogonality relation for these polynomials. In order to interpret the second relation, 
we consider the polynomials Pjn{xn) {{xn} is the set of points, on which the polynomials are 
orthogonal) as functions of m. In this way one obtains one or two sets of orthogonal functions, 
which are expressed in terms of a dual set of (7-orthogonal polynomials (which corresponds to 
the indeterminate moment problem). The second relation in (1.1) leads to the orthogonality 
relations for these dual g-orthogonal polynomials. 

Since this set of g-orthogonal polynomials corresponds to the indeterminate moment prob- 
lem, there are infinitely many orthogonality relations. Using the pair of operators {A, I) and 
the notion of duality, one is able to find only one orthogonality relation (which is dual to 
the orthogonality relation for the corresponding set of polynomials, associated with the de- 
terminate moment problem). Sometimes a measure, which corresponds to this orthogonality 
relation, is extremal and sometimes it is not extremal. It depends on a concrete pair of dual 
sets of polynomials. 

Throughout the sequel we always assume that g is a fixed positive number such that q < 1. 

We use (without additional explanation) notations of the theory of special functions and the 
standard q-analysis (sec, for example, [GR] and [AAR]). We shall also use the well-known 
shorthand notation (ai, • • • , a^; g)„ := {ai;q)n • • • {ak'-,Q)n- 

2. ORTHOGONALITY MEASURES AND DUALITY 

2.1. Orthogonal polynomials, Jacobi matrices and the moment problem 

Orthogonal polynomials are closely related to operators represented by a Jacobi matrices. 
In what follows we shall use only symmetric Jacobi matrices and the word "symmetric" will 
be often omitted. By a symmetric Jacobi matrix we mean a (finite or infinite) symmetric 
matrix of the form 

/60 ao •••\ 
ao 61 ai 

ai 62 fl2 
02 63 03 



M 



(2.1) 



We assume below that all 7^ 0, i = 0, 1, 2, • • •; then are real. Let L be a closed symmetric 
operator on a Hilbert space H-, representable by a Jacobi matrix M. Then there exists an 
orthonormal basis e„, n = 0, 1, 2, • • •, in H, such that 

where e_i = 0. Let \x) = Yl^=oPn{x)en be an eigenvector^ of L with an eigenvalue x, that 
is, L\x) = x\x). Then 



n=0 



n=0 



Equating coefficients of the vector e„, one comes to a recurrence relation for the coefficients 

Pnix): 

anPn+l{x) + bnPn{x) + an-lPn-l{x) = XPn{x). (2.2) 

^Observe that eigenvectors of L may belong to either the Hilbert space Ti or to some extension of Ti.. (For 

example, ii Ti. — LF'{—oq, 00) and in place of L we have the operator d/dx, then the functions e'^'', which do 
not belong to L^{— 00,00), are eigenfunctions of d/dx.) Below we act freely with eigenvectors, which do not 
belong to Ti., but this can be easily made mathematically strict. 
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Since = 0, by setting po{x) = 1 we see that pi{x) = aox — bo/aQ. Similarly we can find 

uniquely Pnix), n = 2, 3, • • •. Thus, the relation (2.2) completely determines the coefficients 
Pn{x)- Moreover, the recursive computation of Pnix) shows that these coefficients Pn{x) are 
polynomials in x of degree n. Since the coefficients a„ and 6„ arc real (because the matrix M 
is symmetric), all coefficients of the polynomials Pn{x) themselves are real. 

Well-known Favard's characterization theorem for polynomials Pn{x), n = 0, 1,2, of 
degree n states that if these polynomials satisfy a recurrence relation 

AnPn+l{x) + BnPn{x) + CnPn-l{x) = xPn{x) 

and the conditions AnCn+i > are satisfied, then these polynomials are orthogonal with 
respect to some positive measure. It is clear that the conditions of Favard's theorem are 
satisfied for the polynomials Pn{x) because in this case the requirements simply reduce to 
inequalities > for n = 0, 1,2, This means that the polynomials Pn{x) from (2.2) 
are orthogonal with respect to some positive measure Ijl{x). It is known that orthogonal 
polynomials admit orthogonality with respect to either unique positive measure or with respect 
to infinitely many positive measures. 

The polynomials Pn{x) are very important for studying properties of the closed symmetric 
operator L. Namely, the following statements are true (see, for example, [Bcr] and [Sim]): 

I. Let the polynomials Pn{x) are orthogonal with respect to a unique orthogonality measure 

j Pmix)pnix)dfJ,{x) = Smn, 

where the integration is performed over some subset (possibly discrete) of M, then the closed 
operator L is self-adjoint. Moreover, the spectrum of the operator L is simple and coincides 

with the set, on which the polynomials Pn{x) are orthogonal (recall that we assume that all 
numbers a„ are non- vanishing) . The measure //(x) determines the spectral measure for the 
operator L (for details see [Ber], Chapter VII). 

II. Let the polynomials Pn{x) are orthogonal with respect to infinitely many different 
orthogonality measures jjl. Then the closed symmetric operator L is not self-adjoint and 
has deficiency indices (1, 1), that is, it has infinitely many (in fact, one-parameter family 
of) self-adjoint extensions. It is known that among orthogonality measures, with respect to 
which the polynomials are orthogonal, there are so-called extremal measures (that is, such 
measures that a set of polynomials {pn{x)} is complete in the Hilbert space with respect 
to the corresponding measure; sec subsection 2.3 below). These measures uniquely determine 
self-adjoint extensions of the symmetric operator L. There exists one-to-one correspondence 
between essentially distinct extremal orthogonality measures and self-adjoint extensions of 
the operator L. The extremal orthogonality measures determine spectra of the corresponding 
self-adjoint extensions. 

The inverse statements are also true: 

I'. Let the operator L be self-adjoint. Then the corresponding polynomials Pnix) are 
orthogonal with respect to a unique orthogonality measure /x, 

j Pm{x)pn{x)dn{x) = Smn, 

where the integral is taken over some subset (possibly discrete) of M, which coincides with the 
spectrum of L. Moreover, a measure /x is uniquely determined by a spectral measure for the 
operator L (for details see [Ber], Chapter VII). 
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II'. Let the closed symmetric operator L be not self-adjoint. Since it is representable 
by a Jacobi matrix (2.1) with 7^ 0, n = 0, 1,2, it admits one-parameter family of 
self-adjoint extensions (see [Ber], Chapter VII). Then the polynomials Pnix) are orthogonal 
with respect to infinitely many orthogonality measures fi. Moreover, spectral measures of 
self-adjoint extensions of L determine extremal orthogonality measures for the polynomials 
{pn{x)} (and a set of polynomials {pn{x)} is complete in the Hilbert spaces L'^ifi) with respect 
to the corresponding extremal measures /x). 

On the other hand, with the orthogonal polynomials Pn{x), n = 0, 1, 2, • • •, the classical 
moment problem is associated (sec [ST] and [Akh]). Namely, with these polynomials (that 
is, with the coefficients and 5^ in the corresponding recurrence relation) real numbers Cji, 
n = 0, 1, 2, • • •, are associated, which determine the corresponding classical moment problem. 
(The numbers Cn are Tiniqucly determined by a,„ and &„.) The definition of the classical 
moment problem consists in the following. Let a set of real numbers Cn, n = 0, 1,2, • • •, be 
given. We are looking for a positive measure /Lt(a;), such that 



where the integration is taken over M. (In this case we deal with the Hamburger moment 
problem.) There are two principal questions in the theory of moment problem: 

(i) Does there exist a measure iJ,{x), such that relations (2.3) are satisfied? 

(ii) If such a measure exists, is it determined uniquely? 

The answer to the first question is positive, if the numbers c^, n = 0, 1, 2, • • •, are those, 
which correspond to a family of orthogonal polynomials. Moreover, a measure /i(a;) then 
coincides with the measure, with respect to which these polynomials are orthogonal. 

If a measure fj, in (2.3) is determined uniquely, then we say that we deal with the deter- 
minate moment problem. In particular, it is the case when the measure fi is supported on a 
bounded set. If a measure, with respect to which relations (2.3) hold, is not unique, then we 
say that we deal with the indeterminate moment problem. In this case there exist infinitely 
many measures ii{x) for which (2.3) take place. Then the corresponding polynomials are 
orthogonal with respect to all these measures and the corresponding symmetric operator L is 
not self-adjoint. In this case the set of solutions of the moment problem for the numbers {c„} 
coincides with the set of orthogonality measures for the corresponding polynomials {pn{x)}. 

Observe that not each set of real numbers Cn, n = 0, 1,2, • • •, is associated with a set of 
orthogonal polynomials. In other words, there are sets of real numbers c„, n = 0, 1, 2, • • • , such 
that the corresponding moment problem does not have a solution, that is, there is no positive 
measure ^u, for which the relations (2.3) arc true. But if for some set of real numbers c^^, 
n = 0, 1, 2, • • •, the moment problem (2.3) has a solution fi, then this set corresponds to some 
set of polynomials Pnix), n = 0, 1, 2, • • •, which are orthogonal with respect to this measure 
fi. There exist criteria indicating when for a given set of real numbers Cn, n = 0,1,2,---, 
the moment problem (2.3) has a solution (see, for example, [ST]). Moreover, there exist 
procedures, which associate a collection of orthogonal polynomials to a set of real numbers 
c„, n = 0, 1, 2, • • •, for which the moment problem (2.3) has a solution (see, [ST]). 

Thus, we see that the following three theories are closely related: 

(i) the theory of symmetric operators L, representable by a Jacobi matrix; 

(ii) the theory of orthogonal polynomials in one variable; 

(iii) the theory of classical moment problem. 





(2.3) 
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2.2. Extremal orthogonality measures 



To a set of orthogonal polynomials Pn{x), n = 0, 1, 2, • • •, associated with an indeterminate 
moment problem (2.3), there correspond four entire functions A{z), B{z), C{z), D(z), which 
are related to appropriate orthogonality measures /x for the polynomials by the formula 

^ A{z) - a{z)C{z) ^ r <W)_ .24^ 
^ B{z)-a{z)D{z) J^^z-t 

(see, for example, [ST]), where a{z) is any analytic function. Moreover, to each analytic 
function a{z) (including cases of constant a{z) and a{z) = ±oo) there corresponds a single 
orthogonality measure /i(t) = Ha{i) and, conversely, to each orthogonality measure fi there 
corresponds an analytic function a such that formula (2.4) holds. There exists the Stieltjes 
inversion formula, which converts the formula (3.1). It has the form 

[fiih + 0) + ii{h - 0)] - [/i(to + 0) + n{to - 0)] 

= lim f\F{t + is)-F{t-ie)]dt). (2.5) 

Thus, orthogonality measures for a given set of polynomials Pn{x), n = 0, 1,2,---, in 
principle, can be found. However, it is very difficult to evaluate the functions A{z), B(z), 
C{z), D{z). In [IM] they are evaluated for particular example of polynomials, namely, for the 

-continuous Hcrmite polynomials hn{x\q). So, as a rule, for the derivation of orthogonality 
measures other methods are used. 

The measures iiaif)-, corresponding to constants a (including a = ±oo), are called extremal 
measures (some authors, following the book [Akh] , call these measures A?^-extremal) . All other 
orthogonality measures arc not extremal. 

The importance of extremal measures is explained by Riesz's theorem. Let us suppose 
that a set of polynomials Pn{x), n = 0, 1,2, • • •, associated with the indeterminate moment 
problem, is orthogonal with respect to a positive measure n (that is, n is a solution of the 
moment problem (2.3)). Let L'^{fJ,) be the Hilbert space of square intcgrablc functions with 
respect to the measure /n. Evidently, the polynomials Pn{x) belong to the space L'^{ijl). Riesz's 
theorem states the following: 

Riesz's theorem. The set of polynomials Pn{x), n = 0, 1, 2, • • is complete in the Hilbert 
space L'^{ij) (that is, they form a basis in this Hilbert space) if and only if the measure fj, is 
extremal. 

Note that if a set of polynomials Pn{x), n = 0, 1,2,---, corresponds to a determinate 
moment problem and ^ is an orthogonality measure for them, then this set of polynomials is 
also complete in the Hilbert space L'^{n)- 

In particular, Riesz's theorem is often used in order to determine whether a certain or- 
thogonality measure is extremal or not. Namely, if we know that a given set of orthogonal 
polynomials, corresponding to an indeterminate moment problem, is not complete in the 
Hilbert space L^(//), where is an orthogonality measure, then this measure is not extremal. 

Note that for applications in physics and in functional analysis it is of interest to have 
extremal orthogonality measures. If an orthogonality measure ^ is not extremal, then it is 
important to find a system of orthogonal functions {fmix)}, which together with a given set 
of polynomials constitute a complete set of orthogonal functions (that is, a basis in the Hilbert 
space LP'{ijl)). Sometimes, it is possible to find such systems of functions (see, for example, 
[CKK]). 
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Extremal orthogonality measures have many interesting properties [ST]: 

(a) If iJ,a{x) is an extremal measure, associated (according to formula (2.4)) with a number 
(T, then iia{x) is a step function. Its spectrum (that is, the set on which the corresponding 
polynomials Pn{x), n = 0, 1,2, are orthogonal) coincides with the set of zeros of the 
denominator B{z) — aD{z) in (2.4). The mass, concentrated at a spectral point xj (that is, 
a jump of ^aix) at the point ,Tj), is equal to (^J^q \Pn{x ■ 

(b) Spectra of extremal measures are real and simple. This means that the corresponding 
self-adjoint operators, which are self-adjoint extensions of the operator L, have simple spectra, 
that is, all spectral points are of multiplicity 1. 

(c) Spectral points of two different extremal measures fJ,a{x) and jia'ix) are mutually 
separated. 

(d) For a given real number xq, always exists a (unique) real value cr, such that the measure 
Haix) has xo as its spectral point. The points of the spectrum of fJ-a{x) are analytic monotonic 
functions of a. 

It is difficult to find all extremal orthogonality measures for a given set of orthogonal 
polynomials (that is, self-adjoint extensions of a corresponding closed symmetric operator). 
As far as we know, at the present time they are known only for one family of polynomials, 
which correspond to indeterminate moment problem. They are the g~ ^-continuous Hermite 
polynomials hn{x\q) (see [IM]). 

If extremal measures /la are known then by multiplying by a suitable factor (depend- 
ing on a) and integrating it with respect to cr, one can obtain infinitely many continuous 
orthogonality measures (which are not extremal). 

2.3. Dual sets of orthogonal polynomials 

A notion of duality for two families of polynomials, orthogonal on finite sets of points, is 
well known. Namely, let Pn{x), n = 0,1,2, - ■■ ,N, he orthogonal polynomials with orthogo- 
nality relation 

TV TV 

(2.6) 

m=0 s=0 

where n,n' = 0,1,2, N, Ws > is a jump of the orthogonality measure in the point Xs, 
and 

n 
k=0 

(ofc and Cn are coefficients in the three-term recurrence relation xpn{x) = anPn+i{x)+hnPn{x)+ 
CnPn-i{x) for the polynomials Pn{x)). Then the dual orthogonality relation is of the form 

JV TV 
'^Pn{Xm)Pn{Xm')v 

E^- (2.7) 

n=0 s=0 

where m,m' = 0,1,2, N (see, for example, [GR], Chapter 7). If one considers the Pn{xm) 
as functions of n, in many cases these functions turn out to be either polynomials in n or poly- 
nomials in u{n), where z/(n) is some function of n. Then the polynomials Pm{n) := Pn{xm), 
m = 0,1,2, - ■ ■ ,N (respectively, Pm{i^{n)) := Pnixm), = 0,1,2, - - - ,N) are orthogonal 
polynomials of n (respectively of iy{n)), for which (2.7) is an orthogonality relation. The 
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polynomials Pj„(i^(™))) = 0, 1, 2, • • • , A'^, are called dual polynomials with respect to the 
Pn{xm), n = 0,1,2, - ■■ ,N. If the dual polynomials {Pm{i^{n))} coincide with the {^'n(w^)}, 
then the polynomials {Pn(7Ti)} are called self-dual. For instance, Racah polynomials and 
g-Racah polynomials both represent families of self-dual polynomials. 

It is not obvious how to extend the notion of duality to polynomials, orthogonal on count- 
able sets of points. In the case of polynomials, orthogonal on a finite set of points, the 
orthogonality (2.7) readily follows from the orthogonality (2.6). Namely, the orthogonality 
(2.6) means that the real (A'^ -|- 1) x (A'' -|- 1) matrix (amn)m,n=o ^^^^ matrix elements 

- orthogonal. Orthogonality of its columns is equivalent to the relation (2.6). Orthogonality 
by rows for the matrix {amn)m,n=o yields the relation (2.7). 

In the case, when we have orthogonality of polynomials on a countable set of points a 
similar conclusion can be false. Let Pn{x), n = 0, 1, 2, • • •, be a set of orthogonal polynomials 
with orthogonality relation 

oo 

XI Pn{Xm)Pn'{Xm)Wm = K^nn' , (2-8) 
m=0 

where n,n' = 0, 1,2,--- and /i„ are some constants. Again, one may consider Pn{xm) as 
functions of n. We are interested in the cases when these functions are polynomials either in 
n or in some v{n). So the question arises: When the dual relation to (2.8), namely, 

oo 

^Pn{Xm)Pn{Xm'))K^ =W~^6mm', (2-9) 

n=0 

is an orthogonality relation for the dual polynomials Pm(z/(n)) := Pn{xrn), m = 0, 1, 2, • • •? It 
follows from Riesz's theorem that this is the case, when the orthogonality measure in (2.8) 
corresponds to determinate moment problem or when this measure corresponds to indetermi- 
nate moment problem and it is extremal. Namely, in both these cases the matrix (amn)m,n=o 
with amn = {h~^WmY^'^Pn{xm) is Orthogonal, that is, 

m 

It is natural to call the polynomials Pm{i^{n)) dual to the polynomials Pn{m). The orthogo- 
nality relation for them is 

oo 
n=Q 

However, very often a function z^(n), such that Pm{i^{n)) := Pnixm), ™ = 0, 1, 2, • • •, are 
polynomials in u{n), does not exist. Nevertheless, sometimes it turns out that there are 
some m-independent bn, n = 0, 1,2, such that Pm{i^{n)) := bnPn{xm), = 0)l5 2, 
are polynomials in iy{n) for an appropriate function iy{n). When the orthogonality measure 
in (2.8) corresponds to determinate moment problem or when this measure corresponds to 



n 
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indeterminate moment problem and it is extremal, then we have the orthogonality relation 
(2.9) for the functions Pm(z/(n)) ■=Pn{xm), which is equivalent to the orthogonality relation 

oo 

^ PmHn)) Pm'{iy{n))b-^h-' = W-^Smm' (2-10) 
n=0 

for the polynomials Pm{i^{n))- In this case it is also natural to call the polynomials Pm{i^{n)) = 
bnPnixm) dual to the orthogonal polynomials pn{m). 

The situation can be sometimes more complicated. Namely, the orthogonality relation for 
polynomials with a weight function, supported on a countable set of points, may be of the 
form (for instance, for the big g-Jacobi polynomials) 

oo oo 

Pn{Xm)Pn'{Xm)Wm + ^ Pn{ym)Pn' {ym)w'^ = hnKn' ■ (2.11) 
m=0 m=0 

Let hn, 6^, i>(n) and v'in) be such functions of n that Pm{i^{Ti)) : — &nPn(^m) Sind P^iy' {n)) : — 
h'^PniUm) are polynomials in v{n) and v'{n), respectively. When the orthogonality measure 
in (2.11) corresponds to determinate moment problem or when this measure corresponds to 

f {a )°° _ \ 

indeterminate moment problem and it is extremal, then the matrix ^"^m,n=o ^ with two 



in' 

\mn/m,n— 

infinite matrices (placed one over another) with matrix elements amn = {h~^WnY^'^Pn{xm) 
and = {h~^w'^^/'^Pn{ym)i is orthogonal, that is, 

m m 
n n n 

Orthogonality of columns of this matrix gives the orthogonality relation (2.11). The orthog- 
onality of rows gives orthogonality of the polynomials Pmiy\ny) and P^{iy'{n)): 



J2 Pm{l^{n))Pm'{y{n))h-^h-^ = W^5mm' (2-12) 

n=0 

oo 

PU^'{n))PU^'{n))h'-\-^ = w'-^b^^,, (2.13) 



n=0 

oo 



Y,P^{v{n))P'^\v\n))h-^^\J = 0. (2.14) 

n=0 

In this case both sets of the polynomials Pmii^'in)) and regarded on a par as 

duals to the set of orthogonal polynomials Pn{x). (For the big g-Jacobi polynomials, these two 
dual sets turn out to be polynomials of the same type, but with different values of parameters; 
see section 4.) If the orthogonality relation (2.11) would contain r terms, then we had r dual 
sets of functions. 

Thus, if we have a set of polynomials {pnix)}, orthogonal on a countable set of points, and 
they correspond to a determinate moment problem or to an indeterminate moment problem 
and the orthogonality measure is extremal, then one can find the corresponding orthogonality 
measure for dual set of polynomials, if they exist. 
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The main goal of this review is to discuss a method of constructing orthogonahty measures 
for a given set of polynomials and their duals in a straightforward manner. This method is 
based on the use of two closed symmetric (self-adjoint) operators, representable (in some 
bases) by Jacobi matrices. In the following sections we shall illustrate how this method works 
by considering families of g-orthogonal polynomials from the Askey scheme. 

Let us emphasize that there are already known theorems on dual orthogonality properties 
of polynomials, whose weight functions are supported on an infinite set of discrete points (see, 
for example, [E], [Sz], [BC] and [Ism]). But it is essential that in most cases (especially in 
the cases of g-polynomials) dual objects are represented by orthogonal functions. Therefore, 
one still needs to make one step further in order to single out an appropriate family of dual 
polynomials from these functions (in those cases when it turns out to be possible). We 
show that this step can be made by choosing the 6„ for the dual polynomials Pm{i^{n)) : = 
bnPn{xm)- Besides, when one considers some dual set with respect to a given family of 
orthogonal polynomials, it is also necessary to investigate the problem of completeness for 
this dual object. In our approach, based on the use of two particular operators, the problem 
of completeness is resolved automatically. 

2.4. List of dual sets of g-orthogonal polynomials 

In this subsection we give a list of dual sets of g-polynomials, orthogonal on countable sets 
of points. Each of these dual pairs will be considered in detail in the subsequent sections. In 
particular, orthogonality relations for them will be explicitly derived. 



g-polynomials 
(determinate moment problem) 


their duals 
(indeterminate moment problem) 


little g-Jacobi 


dual little g-Jacobi 


big g-Jacobi 


dual big g-Jacobi (two sets) 


discrete (?-ultraspherical 


dual discrete g-ultraspherical 


big Q'-Laguerre 


g-Mcixner (two sets) 


alternative g-Charlier 


dual alternative g-Charlier 


Al-Salam-Carlitz I 


g-Charlier (two sets) 


little g-Laguerre 


Al-Salam-Carlitz II 



Let us exhibit these dual pairs explicitly. 

Little g-Jacobi polynomials and their duals. Little g-Jacobi polynomials, given by the 
formula 

Pn{X;a,b\q) := 20i(g'",a&g"+^ aq; q,qX), (2.15) 

are orthogonal for < a < q~^ and b < q~^. The dual little g-Jacobi polynomials, corre- 
sponding to the polynomials (2.15) with the same values of the parameters a and b, are given 
as 

d„(/x(m);a,6|g) := 3<^i(g-"*, a6g™+\ g""; bq; q,q-/a), (2.16) 

where ju(m) = g"™ + abq"^^^. Since these polynomials are absent in the Askey g-scheme 
[KSw], we give the orthogonality relation for these polynomials: 



m=0 



^ [1- abq){aq,q;q)r- 
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The polynomials dn{fi{m)) correspond to the indeterminate moment problem and the ortho- 
gonality mcasTirc here is extremal. The duality of these polynomials to the set of the little 
g'-Jacobi polynomials was first observed in [AKl] (see also [AK3]). 

Big g-Jacobi polynomials and their duals. Big g-Jacobi polynomials, given by the 
formula 

Pn{X;a,b,c;q) := 3^2(9"", a6g"+\ A; aq,cq; q,q), (2.17) 

are orthogonal for < a, 6 < q~^ and c > 0. The dual big g-Jacobi polynomials, associated 
with the polynomials (2.17) with the same values of the parameters a,b,c, are given as 

Dniiiim); a, b, c\q) := sMq'"", abq^+\q-''; aq, abq/c; q, aq^'+'/c), (2.18) 

where /i(m) = q~"^ + abq"^~^^. The second set of dual polynomials with respect to (2.17) is 
obtained from the polynomials (2.18) by the replacements a,b,c — > b,a,ab/c, respectively. 
Again, since these polynomials are absent in the g-Askey scheme, we give here the orthogo- 
nality relation for the polynomials (2.18): 



E°° {1 -abq^"^+'^){aq,abq, abq/c; q)m f , .m m(m-i)/2 n ft \\ n t t \\ 
{abq'^, c/ a; q) 00 {aq/c,q;q)n 



m=0 (^bq){bq,cq,q;q) 



{bq, cq;q)oo {aq, abq/ c; g)„g'* 

The polynomials Dn{iJ,{m)) correspond to the indeterminate moment problem and the or- 
thogonality measure here is not extremal. 

Discrete g-ultraspherical polynomials and their duals. Discrete g-ultraspherical poly- 
nomials C^\x;q), a > 0, are a particular case of the big g-Jacobi polynomials 

Ci''^\x;q) = Pnix : a,a,~a;q) = 3(f)2iq''\a'^q"'^^,x; aq,-aq; q,q). (2.19) 

An orthogonality relation for Cn\x; q) follows from that for the big g-Jacobi polynomials and 
it holds for positive values of a. We can consider the polynomials Cn\x\q) also for other 
values of a. In particular, they are orthogonal for imaginary values of a and x. In order to 
dispense with imaginary numbers in this case, the following notation is introduced: 

Cf\x;q) := (-i)"C(-"')(ix; g) = (-i)"P„(ix; ia, ia, -ia; g), (2.20) 

The orthogonality relation for them is of the form 

— / o.„2^ <-2fc+i(.va5 ,q)L^^,^^Waq ,q} 

s=0 W .9 ;5 

^ (-gg^ q^)oo (1 + aq) a^^+^ (g; g)2fc+i ^(fc+2)(2fc+i) 
(9; 9^)00 (l-|-ag4fe+3) {-aq;q)2k+i 



The formula 

/ ^—x „2„x+l „—n 

g,-g"+M, (2.21) 



D':^\li{x\a^)\q) := Dn{n{x;a^);a,a,-a\q) := ^(t)2' 



aq, —aq 
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where a^) = + a^g^"'"^ and Dn{iJ,{x; a^)) are dual big g-Jacobi polynomials, gives dual 
discrete g-ultraspherical polynomials. They correspond to indeterminate moment problem. 
The dual orthogonality relation for them (when > 0) follows from the orthogonality relation 
for dual big g-Jacobi polynomials. 

For the polynomials \iJ,{x;a^)\q) with imaginary a we have 

-D^"^)(/x(x;-a^)|g) := DnifJ-{x;-a^);ia,ia,-ia\q) 

q,-q^+'). (2.22) 



3^2 



g ^,-a^g^+-^,g " 



lag, —lag 



~ (a^) 

These polynomials are dual to the polynomials Ck {x; q) from (2.20). In this case there are 
also infinitely many orthogonality relations, which are considered in section 5. 

Big g-Laguerre polynomials and g-Meixner polynomials. Big g-Laguerre polynomials, 
given by the formula 

P„(A;a,6;g) := 3^2(g"", 0, A; ag,6g; g,g), 

are orthogonal for < a < q~^ and 6 < 0. The dual polynomials coincide with g-Meixner 
polynomials Mn{q~^; a, —b/a; q) and Mn{q~^; b, —a/b; g), where 

M„(g-^;a,6;g) := 2<^i(g~'', g""^; ag; q,-q''+^/b). 

We obtain orthogonality relations for the g-Meixner polynomials Mn{q~^;a,b;q) with b < 
and 5 > in section 6. 

The duality relation between big g-Laguerre polynomials and g-Meixner polynomials was 
studied in [AAK]. The appearance of g-Meixner polynomials as a dual family with respect to 
the big g-Laguerre polynomials is quite natural because the transformation g q~^ interre- 
lates these two sets of polynomials, that is, 

Mn{x; b, c; g"^) = (g^"/6; g)„ Pn{qx/b; 1/5, -c; g). 

Alternative g-Charlier polynomials and their duals. Alternative g-Charlier polynomials 
are given by the formula 

Kn{X; a; g) := 2<^i(g-'^, -ag"; 0; g, gA) . (2.23) 
They are orthogonal for a > 0. Their duals are the polynomials 

d„(Mm);a;g) := 3</>o(g-™, -og"*, g""; -; g, -g7«) , M"^) :=g-™-ag-, (2.24) 

which correspond to the indeterminate moment problem (see [AK5]). They are also absent in 
the g-Askey scheme. The orthogonality relation for these polynomials is 



^0 ("'^^ ; 9)00 (g; q)m a"g"i"+i)/^ 



Al-Salam— Carlitz I polynomials and g-Charlier polynomials. Al-Salam-Carlitz I poly- 
nomials, given by the formula 

Ui-\x- q) :=(-a)V(«-«/S0i(g-",x-^ 0; g;xg/a). 
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are orthogonal for a < 0. There are two sets of dual polynomials [KK]. They coincide with 
two sets of g-Charlier polynomials Cn{q~^', —a; q) and Cn{q~^', — q), where 

Cn{q-''; a; q) := 2</>i(g-", g'^; 0; g;-(z"+V«)- 

Little g-Laguerre polynomials and Al-Salam— Carlitz II polynomials. The little q- 
Laguerre polynomials are given by the formula 

Pn{x; a\q) := 2<?^i(g~",0; aq; gx) = (a~^g~"; g)~S'?!'o(g~", -; q;x/a). 

They are orthogonal for < a < The dual polynomials with respect to them are the 

Al-Salam-Carlitz II polynomials (see [AK2]) 

V^''\x;q) ■.= {-aTq~<'^-^y\UQ-\x- -; q-q^/a). 



3. LITTLE g-JACOBI POLYNOMIALS AND THEIR DUALS 
3.1. Pair of operators {Ii,J) 

Let K be a separable complex Hilbert space with an orthonormal basis /„, n = 0, 1, 2, ■ • •. 

The basis determines uniquely a scalar product in 7Y. In order to deal with a Hilbert space of 
functions on a real line, we fix a real number a such that < a < q~^ and realize our Hilbert 
space in such a way that basis elements are monomials: 

fn = fnix) CnX , 

where 

ci = l, ci = g(^-^'W^^^^, n= 1,2,3,.... 

{q; q)n 

Thus, in fact, our Hilbert space depends on the number a and can be denoted as TYa- 

We fix two real parameters a and b such that b < q~^, < a < q~^, and define onH = Ha 
two operators. The first one, denoted as q'^° and taken from the theory of representations of 
quantum group [/g(sui,i), acts on the basis elements as 

q''fn = iqa)'/^qVn. (3.1) 
The second operator, denoted as Ii, is given by the formula 

h fn = —O'nfn+l — ^n-l/n-l + ^n/n; (3-2) 



where 



^ ^ ^1/2 n+1/2 - - Qg"+^)(1 - 6g»+^)(l - ab^) 

^ (1 - a6g2n+2)y(l _ a692n+l)(l _ a6g2n+3) 

/(I - ag"+i)(l - a?)g"+M (1 - g")(l - 6g") 

+ a, 



1 - abq^'^+^ \ 1 - a6g2«+2 ' i _ a6g2n 

The expressions for a„ and 6„ are well defined. The operator Ii is symmetric. 

Since a„ ^ and 6„ — > when n oo, the operator Ii is bounded. Therefore, we assume 
that it is defined on the whole space H. For this reason, I\ is a self-adjoint operator. Let us 
show that I\ is a Hilbert-Schmidt operator (we remind that a bounded self-adjoint operator 
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is a Hilbert-Schmidt operator if a sum of its squared matrix elements in an orthonormal basis 
is finite; the spectrum of such an operator is discrete, with a single accumulation point at 0). 
For the coefficients a„ and 6„ from (3.2), we have 

Qn+i/on q, bn+i/bn ^ Q when n ^ oo. 

Therefore, for the sum of all matrix elements of the operator Ii in the canonical basis we have 
^^(2a„ + 6„) < oo. This means that /i is a Hilbert-Schmidt operator. Thus, the spectrum of 
Ii is discrete and has a single accumulation point at 0. Moreover, a spectrum of Ii is simple, 
since Ji is representable by a Jacobi matrix with a„ 7^ (see [Ber], Chapter VII). 
To find eigenfunctions S,x{x) of the operator Ii, Ii^x{x) = ^^x{x), we set 

00 

n=0 

Acting by the operator Ii upon both sides of this relation, one derives that 

00 oo 

/3n(A) {anfn+1 + On-l/n-l - bnfn) = PnWfn, 
n=0 n=0 

where a„ and 6„ are the same as in (3.2). Collecting in this identity all factors, which multiply 
/„ with fixed n, one derives the recurrence relation for the coefficients PnW- 

/3„+i(A)a„ + /3„_i(A)a„_i - /?„(A)6„ = -A/3„(A). 

The substitution 

reduces this relation to the following one 

An/3'n+lW + CnP'n-l{X) - (An + Cn)P'nW = -XP'nW 

with 

_ g"(l-ag"+i)(l-a6g"+^) _ ag"(l-g")(l-6g") 

"~ (l-a6g2n+l)(l_a6g2n+2)' - (l_a;,^2n)(l_(,5^2n+l)- 

This is the recurrence relation for the little g-Jacobi polynomials 

Pn{X;a,b\q) := 20i(g~"',a&g"'^^ aq; q,qX) (3.3) 
(see, for example, formula (7.3.1) in [GR]). Therefore, /3^(A) = Pn{X;a,b\q) and 

^2n+l^\ 1/2 



/9n(A) 



f {abq,aq;q)n{l-abq^''+^) \ 
-77 r-r. , w pn{\;a,b\q). (3.4) 

V (&g,g;g)n(i-aMH) / 



For the eigenfunctions ^\{x) we have the expression 

/(a6g,ag;g)„ (1 - a&g2"+i)\ V2 , , x . / x 
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Since the spectrum of the operator Ii is discrete, only a discrete set of these functions belongs 
to the Hilbert space H. This discrete set of functions determines a spectrum of Ii. 

Now we look for a spectrum of the operator Ii and for a set of polynomials, dual to the 
little g-Jacobi polynomials. To this end we use the action of the operator 

J := (ga)V2g--^o + {qa)-^/'^abq'^°+^ 

upon the eigenfunctions ^a(^)) which belong to the Hilbert space H. In order to find how this 
operator acts upon these functions, one can use the g-difference equation 

(g-'^ + a6g"+i)p„(A) = aA-^fegA - l)p„(gA) + X'^l + a)p„(A) + X'^X - l)pn{q-'X) (3.6) 

for the little g'-Jacobi polynomials Pn{X) = PniX; a,b\q) (see, for example, formula (3.12.5) 
in [KSw]). Multiply both sides of (3.6) by dn fn{x) and sum up over n, where d„ are the 
coefficients of Pn{X;a,h\q) in the expression (3.4) for the l3n{X). Taking into account the first 
line in formula (3.5) and the fact that Jfn{x) = (g"" + a6g"+^)/„(a;), one obtains the relation 

JUx) = aX-\bqX - 1) ^,a(x) + A-i(l + a) U^) + A'^A - 1) Vul^^)- (3.7) 

It will be shown in the next section that the spectrum of the operator Ii consists of the 
points A = g", n = 0, 1, 2, • • •. Thus, we see that the pair of the operators Ii and J form a 
Leonard pair (see [Terl], where P. Terwilliger has actually introduced this notion in an effort 
to interpret the results of D. Leonard [Leo]; sec also [Ter3], which contains a review on how 
one can employ Leonard pairs to describe properties of orthogonal polynomials). We remind 
to the reader that a pair of operators Ri and i?2, acting on a linear space £, is a Leonard pair 
if 

(a) there exists a basis in C, with respect to which the operator Ri is diagonal, and the 
operator R2 has the form of a Jacobi matrix; 

(b) there exists another basis of JC, with respect to which the operator R2 is diagonal, and 
the operator Ri has the form of a Jacobi matrix. 

Properties of Leonard pairs of operators in finite dimensional spaces arc studied in detail. 
Leonard pairs in infinite dimensional spaces arc more complicated and only some isolated 
results are known in this case (see, for example, [Tcr2]). 

3.2. Spectrum of Ii and orthogonality of little g-Jacobi polynomials 

The aim of this section is to find, by using the Leonard pair (/i, J), a basis in the Hilbert 
space H, which consists of eigenfunctions of the operator Ii in a normalized form, and to 
derive explicitly the unitary matrix U, connecting this basis with the canonical basis /„, 
n = 0, 1,2, in 7^. This matrix directly leads the orthogonality relation for the little 
g-Jacobi polynomials. 

Let us analyze a form of spectrum of the operator Ii. If A is a spectral point of the 
operator /i, then (as it is easy to see from (3.7)) a successive action by the operator J upon 
the function (eigenfunction of Ii) leads to the functions 

Cgmx, m = 0,±l,±2,---, (3.8) 

which are eigenfunctions of Ii with eigenvalues g™A. However, since Ii is a trace class operator, 
not all these points can belong to the spectrum of Ii, since g~™A — > 00 when m ^ 00 if A / 0. 
This means that under a successive action by Ji upon on some step the last term in (3.7) 
must vanish. Thus, under the action by Ii upon for some A' the coefHcient A'— 1 of ^g-i\i {x) 
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in (3.7) vanishes. Clearly, it vanishes when A' = 1. Moreover, this is the only possibility for 
the coefficient of ^q-ix'{x) in (3.7) to vanish, that is, the point A = 1 is a spectral point for 
the operator Ii. Let us show that the corresponding eigenfunction ^i{x) = ^qo{x) belongs to 
the Hilbert space 7i. 

Observe that by formula (II. 6) of Appendix II in [GR], one has 

p„(l;a,6|g) =2<^i(g-",a6g"+i; aq; q,q) = (a6g"+^)". 

[aq; Q)n 

Since (6~^g~";g)n = (6g; g)n(-&~^9~^)"g""^""^^/^ this means that 

Pnil;a,b\q) = f^{-arq-^-+'y'. 
Therefore, due to (3.5) for the scalar product (^i(x),^i(x)) we have 

/c r \ i t w (abq, aQ; q)n (1 - a&g^"+^) 2/. x 

^ [bq, q; g)n (1 - abq) (aq)'' 



n=Q 



^ jabq, bg;g)n(l-a6g^"+^) ^„^„2 ^ (abg^; q)^ 
^ {aq,q;q)n{i-abq) {aq;q)oo 

The last leg of this equality is obtained from formula (A.l) of Appendix. Thus, the series (3.9) 

converges and, therefore, the point A = 1 actually belongs to the spectrum of the operator Ii. 

Let us find other spectral points of the operator Ii (recall that a spectrum of Ii is discrete). 
Setting A = 1 in (3.7), we see that the operator J transforms £,qo{x) into a linear combination 
of the functions ^q{x) and ^qQ{x). Moreover, ^q{x) belongs to the Hilbert space H, since the 
series 

ibq,q;q)n (1 - abq){aq)'^ 

is majorized by the corresponding series for ^qVi{x), considered above. Therefore, ^q{x) belongs 
to the Hilbert space Ti and the point q is an eigenvalue of the operator Ii . Similarly, setting A = 
q in (3.7), we find that {x) is an eigenfunction of Ii and the point belongs to the spectrum 
of Ii. Repeating this procedure, we find that ^q^{x), n = 0, 1, 2, • • •, arc eigenfunctions of Ii 
and the set q^, n = 0, 1, 2, • • •, belongs to the spectrum of /i. So far, we do not know yet 
whether other spectral points exist or not. 

The functions iqn{x), n = 0, 1,2, are linearly independent elements of the space H 
(since they correspond to different eigenvalues of the sclf-adjoint operator Ii). Suppose that 
values g", n = 0, 1, 2, • • •, constitute a whole spectrum of the operator Ii. Then the set of 
functions Cq"(^)) = 0, 1, 2, • • •, is a basis in the Hilbert space TL. Introducing the notation 
S„ := ^qn{x), n = 0, 1, 2, • • •, we find from (3.7) that 

JH„ = -a(?-"(l - 6g"+i) + g-"(a + 1) S„ - ^""(1 - q^) H„_i. (3.10) 

As wc sec, the matrix of the operator J in the basis E!„, n = 0, 1, 2, • • •, is not symmetric, 
although in the initial basis /„, n = 0, 1,2, it was symmetric. The reason is that the 
matrix A = {amn) with entries 

amn ■= Pm{q'^), m, n = 0, 1, 2, • • • , 
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where Pmio."') are the coefficients (3.4) in the expansion ^qn{x) = J2m (^m{<f')fm{x), is not 
unitary. (This matrix connects the bases {/„} and It is equivalent to the statement 

that the basis S„ := Cg"(3^)) = 0, 1,2,---, is not normahzed. To normalize it, one has 
to multiply H„ by corresponding numbers c„ (which are not known at this moment). Let 
= n = 0, 1,2,---, be a normalized basis. Then the matrix of the operator J is 

symmetric in this basis. Since J has in the basis the form 

JS„ = -c-jic„ag-"(l - 6g"+i) + q-^{a + 1) H„ - c-^c„g-"(l - g") (3.11) 

then its symmetricity means that 



that is, Cn/cn-i = ^/aq{l — bq"-)/{l — q^). Therefore, 

Cn = ciaqr/^ 



{bq; qfr!"^ 

1/2 ' 



{q; q)n 

where c is a constant. 

Now instead of the expansion (3.5) we have the expansions 

4n(x) = ^n{x) =Y,CnPm{qnfm{x), (3.12) 
m 

which connect two orthonormal bases in the space H. This means that the matrix (a^n), 
m, n = 0, 1, 2, • • •, with entries 

«™» = c„M<f) = c ( (a,r- (a>.,a,;,).,.(l-a6,'-"^-) \ 

V W;5)n {bq,q;q)„,{l- abq) J 

(3.13) 

is unitary, provided that the constant c is appropriately chosen. In order to calculate this 
constant, we use the relation Y11^=q l^mnP = 1 for n = 0. Then the sum in this relation is a 
multiple of the sum in (3.9) and, consequently. 



1/2 

oo 



(gg; g) 

{abq; q)]^ 



Thus the in (3.12) and (3.13) are real and equal to 

^ ^ f (Qg;g)oo {bq]q)n{aqT V''^ 
\{abq;q)oo {q;q)n J 

The matrix (dmn) with entries (3.13) is orthogonal, that is, 

^ ^ ^mn^m'n ~ ^mm' i ^ ^ O'mn&mn' ~ ^nn' ■ (3-14) 

n m 

Substituting into the first sum over n in (3.14) the expressions for amn, we obtain the identity 

E°° {bq;q)n{aqT , n u \ , n u \ 
— — prain ■,a,b\q)pm'{q ■,a,b\q) 

n=0 
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_ (abg^;g)oo {1 - abq)(aq)'^ {bq,q;q)m ^ ^ ^^-^ 
{aq;q)oo {1- abq'^"'+^){abq,aq-q)m """" ' 

which must yield the orthogonality relation for the little g-Jacobi polynomials. An only gap, 
which appears here, is the following. We have assumed that the points g", n = 0, 1,2, 
exhaust the whole spectrum of the operator Ii. Let us show that this is the case. 

Recall that the self-adjoint operator /i is represented by a Jacobi matrix in the basis /„, 
n = 0, 1,2, • • •. According to the theory of operators of such type (see, for example, [Ber], 
Chapter VII; a short explanation is given in section 2), eigenfunctions of Ii are expanded 
into series in the monomials /„, n = 0, 1,2, • • •, with coefficients, which are polynomials in 
A. These polynomials are orthogonal with respect to some positive measure (i//(A) (moreover, 
for self-adjoint operators this measure is unique). The set (a subset of M), on which the 
polynomials arc orthogonal, coincides with the spectrum of the operator under consideration 
and the spectrum is simple. Let us apply these assertions to the operator Ii. 

We have found that the spectrum of /i contains the points g", n = 0, 1,2, If the 
operator Ii had other spectral points x, then on the left-hand side of (3.15) there would 
be other summands fixf, Pmi^k] ^k) Pm'i^ki '^i ^q)i corresponding to these additional points. 
Let us show that these additional summands do not appear. To this end we set m = m' = 
in the relation (3.15) with the additional summands. Since po{x;a,b\q) = 1, we have the 
equality 

{bq;q)n{aqT sr^ _ {abq^;q)oo 

According to the g-binomial theorem (see formula (1.3.2) in [GR]), we have 

{bq;q)n{aqT {abq'^;q)c 



n=0 



/oo 



{q;q)n {aq;q)c 



(3.16) 



Hence, Hx,. = and this means that additional summands do not appear in (3.15) and it 
does represent the orthogonality relation for the little g-Jacobi polynomials. 

By using the operators /i and J, which form a Leonard pair of infinite dimensional sym- 
metric operators, we thus derived the orthogonality relation for little g-Jacobi polynomials. 

The orthogonality relation for the little g-Jacobi polynomials is given by formula (3.15). 
Due to this orthogonality, we arrive at the following statement: The spectrum of the operator 
Ii coincides with the set of points g", n = 0, 1,2, The spectrum is simple and has one 
accumulation point at 0. 

3.3. Dual little g-Jacobi polynomials 

Now we consider the second identity in (3.14), which gives the orthogonality relation for 
the matrix elements a^n, considered as functions of m. Up to multiplicative factors these 
functions coincide with 

Fn{x\a,b\q) := 2<t>i{x,abq/x\ aq; g,g"+^), (3.17) 
considered on the set x G {g~"* | m = 0, 1, 2, • • •}. Consequently, 

. / (a-q; q)oo {bq;q)n (.n-m (abg, aq;q)rn{l- abq'^"'+^) \ ^'"^ „ , m. „ Un\ 

amn — Tl ^ 7 ^ yaq) 77 ^ ^nVl , «, o\q) 

\{abq;q)oo {q;q)n [bq,q;q)m J 
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and the second identity in (3.14) gives the orthogonahty relation for the functions (3.17): 

{l-abq)iaqribq,q;qU "^^^^ ' 
(a6g^;g)oo (g;g)„(ag)~" ^ 

= —, ^ 77 ^ (>nn'- (3.18) 

Kaq\ q)oo [bq; q)n 

The functions Fn{x]a,h\q) can be represented in another form. Indeed, one can use the 
relation (III. 8) of Appendix III in [GR] in order to obtain that 

F„(g— ;a,6|g) = l^I^«-^^(a6g-+i)-3<^i(g-"^, a6g-+\ g-"; bq- q,q^/a) 

[aq; q)m 

= ("^^'"(^^'g)"^ a"^g"^("'+^)/^3'/'i(g~'",Q&g'"+\g~"; bq; q,q''/a). (3.19) 

{aq;q)m 

The basic hyper geometric function in (3.19) is a polynomial of degree n in the variable 
/i(m) := q~"^ + abq^'^^, which represents a g-quadratic lattice; we denote it as 

dMmy,a,b\q) := ^Mq'"" , (^b q"^+\ q'^^ ; bq; q,q''/a). (3.20) 



Then formula (3.18) yields the orthogonality relation 

-abq'^"'+^){abq,bq 
(1 - abq){aq,q;q) 



E°° [1 - abq^"'+'^){abq,bq;q)m ^ / 



m=0 

{abq^;q)oo {q;q)n{aqy" . „t ^ 

= -7 T 77 ^ <^nn' (3-21) 

{aq; q)oo {bq; q)n 

for the polynomials (3.20). We call the polynomials dn{fJ'{'m);a,b\q) dual little q-Jacobi poly- 
nomials. 

Note that these polynomials can be expressed in terms of the Al-Salam-Chihara polyno- 
mials 



^ . ,i N {ab;q)n , fq '^,az,az ^ 
Qn{x;a,b\q) = — — 3(^)2 



with the parameter q > 1. An explicit relation between them is 

qn{n-l)/2 

dn{n{x);P/a, l/afSq \ q) = Qn{an{x)/2; q, /3|g ). 

(-p) (l/«P;9)n 

Ch. Berg and M. E. H. Ismail studied this type of Al-Salam-Chihara polynomials in [BI] and 
derived complex orthogonality measures for them. But [BI] does not contain any discussion 
of the duality of this family of polynomials with respect to little g-Jacobi polynomials. 

Observe that the dual polynomials (3.20) can be also expressed in terms of the little 
g-Jacobi polynomials (3.3): 

d,(M(m);a,6|g) = i^^i^(-a)-™g— (-+i)/2p^(g"; a, 5|g). 
[bq/c; q)m 

A recurrence relation for the polynomials c/„(/x(m); a, 6|g) is derived from formula (3.6). 
It has the form 

{q-^ + afeg'^+i) dMm)) = - ag-"(l - 6g"+i) dn+Mm)) 
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+ + a) dMm)) - q-^{l - q^) dn-Mm)), 

where dn{lJ^{m)) = ci„(//(m); a, Comparing this relation with the recurrence relation 

(3.69) in [All], we see that the polynomials (3.20) are multiple to the polynomials (3.67) 
in [All]. Moreover, if one takes into account this multiplicative factor, the orthogonality 
relation (3.21) for polynomials (3.20) turns into relation (3.82) for the polynomials (3.67) in 
[All], although the derivation of the orthogonality relation in [All] is more complicated than 
our derivation of (3.21). The authors of [All] do not give an explicit form of their polynomials 
in the form similar to (3.20). Concerning the polynomials (3.67) in [All] see also [Gro]. 

Let be the Hilbert space of functions on the set m = 0, 1, 2, • • • with the scalar product 

/i,/2)=> 7z ^^"7 r a q /i(m)/2(m), 

^ (1 - abq) {aq, q; q)r- 



m=0 



where weight function is taken from (3.21). The polynomials (3.20) are in one-to-one corre- 
spondence with the columns of the unitary matrix (a„m) and the orthogonality relation (3.21) 
is equivalent to the orthogonality of these columns. Due to (3.14) the columns of the matrix 
(omn) form an orthonormal basis in the Hilbert space of sequences a = {a„ | n = 0, 1, 2, • • •} 
with the scalar product (a, a') = Yln^na'n- This assertion is equivalent to the following one: 
the set of polynomials dn{ii(m); a, b\q), n = 0, 1, 2, • • •, form an orthogonal basis in the Hilbert 
space l^. This means that the point measure in (3.21) is extremal for the dual little q-Jacobi 
polynomials dn{n{m); a, b\q). 

4. BIG g-JACOBI POLYNOMIALS AND THEIR DUALS 

4.1. Pair of operators (hjJ) 

We fix three real numbers a, b and c such that < a < q~^, < b < q~^, c < and 
consider on the Hilbert space H = Ha, introduced in subsection 3.1, the following symmetric 
operator I2: 

h fn = fln/n+l + An-l/n-l — bnfn, (4.1) 

where 



a _ = f-ac V(l - g")(l - Qg")(l - &g")(l - abq^) (1 - cg")(l - abc-\^ 

_ (l-ag"+^)(l-a6g'^+i)(l-cg"+i) _ (1 - g")(l - 6g")(l - a6g"/c) _ 

"~ (l-a6g2n+i)(i_o6g2n+2) {1 - abq^''){l - abq'^^+^) ' 

This operator is bounded. Therefore, we assume that it is defined on the whole Hilbert space 
H. This means that I2 is a self-adjoint operator. Actually, I2 is a Hilbert-Schmidt operator. 
To show this we note that for the coefficients a„ and 6„ from (4.1) one obtains that 

a„+i/a„ g^/^, bn+i/bn q when n ^ 00. 

Therefore, ^n) < 00 and this means that I2 is a Hilbert-Schmidt operator. Thus, 

the spectrum of I2 is simple (since it is representable by a Jacobi matrix with a„ / 0) , discrete 
and have a single accumulation point at 0. 

To find eigenfunctions tp\{x) of the operator I2, l2tp\{x) = X'tp\{x), we set 

00 

M^) = ^Pn{\)fn{x). (4.2) 
n=0 
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Acting by the operator I2 on both sides of this relation, one derives that 

/3n(A)(a„/n+l + an-lfn-1 - Kfn) = /3n(A)/„, 

n 

where o„ and 6„ are the same as in (4.1). Collecting in this identity factors, which multiply 
with fixed n, we arrive at the recurrence relation for the coefficients /3n(A): 

a-nPn+lW + an-lPn-lW - KPn{>) = A/?n(A). 

Making the substitution 

^n(A) = ( a - \ ^-n(n+3)/4 



Xabq/c, bq,q;q)n (1 - ahq){-ac) 
we reduce this relation to the following one 

^n/3;+i(A) + C„/3;_i(A) - (A„ + Cn- 1K(A) = A/3;(A) 

with 

_ (1 - ag"+^)(l - cg"+^)(l - a6g"+i) 
" ~ (1 - a6g2n+i) (1 _ a6g2n+2) ' 

_ -acg"+i(l - g")(l - 5g")(l - a^c-^g") 
~ (1 - a6g2n) (1 _ a6g2n+i) " 

It is the recurrence relation for the big g-Jacobi polynomials 

Pn(A;a,6,c;g) := 3(^2(9""", afe^'^'^'S A; ag,c5; q,q) (4.3) 

introduced by G. E. Andrews and R. Askey [AA] (see also formula (7.3.10) in [GR]). Therefore, 
l^'ni}) = Pn{>^;a,b,c;q) and 

MX)=( ^ f 7- f - ff"' ) F„(A; ». 6. c; g). (4.4) 

V (a6g/c, 6g, g; g)n (1 - a6g) (-ac)" / 

For the eigenfunctions ipx{x) we have the expansion 

/ / N >fi f {abq, aq, eg; g)„ (1 - abg^^+i) ^^ n . / n 

^ \{abq c,bq,q;q)n (1 - abq)[-ac)'' J 



_ ^ n/4 (ag;g)n ( {abq,cq-q)n{l-abq^'^+^) \ ^'^ n(n+3)/4 p / x . „ z, ^. „^ 
"no (^!^)" V(a^'g/c,bg;g)n(l-a6g)(-ac)"; ^ i^n(A,a,6,c,g)x . 

(4.5) 

Since the spectrum of the operator I2 is discrete, only a discrete set of these functions belongs 
to the Hilbcrt space TL. 

In what follows we intend to study a spectrum of the operator I2 and to find polynomials, 
dual to big g-Jacobi polynomials. It can be done with the aid of the operator 

J := {aqf'^q-^° + {aq)-^/'^abq-^°+\ 
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which has been already used in the previous case in subsection 3.1. In order to determine 
how this operator acts upon the eigenfunctions ipx{x), one can use the g-difFerence equation 

(g-" + a6(?"+i) Pn{X) = aqX~HX - l){bX - c) P„(gA) 

-[X-^acq{l +q)- X-^q{ab + ac + a + c)] P„(A) + X-^{X - aq){X - cq) Pn{q~^X), (4.6) 

for the big g-Jacobi polynomials PniX) = Pn{X;a,b,c;q) (see, for example, formula (3.5.5) in 
[KSw]). Multiply both sides of (4.6) by kn fn{x), where kn are the coefficients of Pn{X; a, b, c; q) 
in the expression (4.4) for the coefficients Pn{X), and sum over n. Taking into account formula 
(4.5) and the fact that J fn{x) = (g"" + afeg""*"^) fn{x), one obtains the relation 

JMx) = aqX-\X -l){bX-c) ^qx{x) 

-[X-'^acq{l + q)- X-'^q{ab + ac + a + c)] i/jxix) + X-'^{X - aq){X - cq) %-ix{x). (4.7) 

It will be shown in the next section that the spectrum of the operator I2 consists of the points 

ag", cq"', n = 0, 1,2, • ■ •. The matrix of the operator J in the basis of eigenfunctions of I2 
consists of two Jacobi matrices (one corresponds to the spectral points ag", n = 0, 1,2, • • •, 
and another to the spectral points eg", n = 0, 1,2, • • •). In this case, the operators I2 and J 
form some generalization of Leonard pair. 

4.2. Spectrum of I2 and orthogonality of big g-Jacobi polynomials 

As in subsection 3.2 one can show that for some value of A (which must belong to the 
spectrum) the last term on the right side of (4.7) has to vanish. There are two such values 
of A: A = aq and A = eg. Let us show that both of these points are spectral points of the 
operator I2. Observe that, according to (4.3), 

P„(ag;a,6,e;g) := 2<^i(g-", a6g"+^ eg; g,g) = (^W; g)" (^^)n gn(n+i)_ 

(eg; g)n 



Therefore, since 
one obtains that 

Likewise, 



(e/a6g";g)„ = (a6g/e; g)„(-e/a6)"g-"("+i)/2 , 

P„(«?; a, h, e; g) := W^(_e)",Mn4-i)/2. (4.3^ 
leg, gjn 

P4eg;a,6,e;g) := fek(-a) ^("+1)/^ 

(ag;g)n 

Hence, for the scalar product {ipaq{x),tpaqix)) we have the expression 

V (1 - (^bq, aq, cq; q)^ n(n+s)/2 p2( ^ ^ ^ . 

il-abq)iabq/c,bq,q;q)n{-ac)n'^ P„(ag,a,^e,gj 

_ (1 - Q&g^''"^^) jabq/c, abq, aq; q)n n(n-i)/2 _ {abq^ ,c/a;q)^ 
^ (1 - a6g)(6g,eg,g;g)„ (-a/e)" {bq,cq;q)oo 

where the relation (A.6) from Appendix has been used. Similarly, for {i/jcq{x) , ipcqix)) one has 
the expression 

V - abq'^'''^^){abq,aq,cq;q)n ^-„(„+3)/2 p^^cq;a,b,c;q) = ("^g^"/^; g)o° ^ (4 iq) 
^ {1 - abq){-ac)'' {abq/c,bq,q;q)n ' ' ' ' {aq, abq /c;q) 00' 
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where formula (A. 7) from Appendix has been used. Thus, the values \ = aq and X = cq are 
spectral points of the operator I2. 

Let us find other spectral points of l2- Setting X = aq in (4.7), we see that the operator 
J transforms ipaqix) into a linear combination of the functions ilj^q2{x) and il^aq{x). We have 
to show that ipaq'^i^) ^-Iso belongs to the Hilbert space 7i, that is, that 



00 



/, / \ {abq,aq,cq;q)n{l - abq^'^+^) -n(n+3)/2 u^f 2 ^ ^/ 

{'4>aq^,i^aq^) = >^ / , / , TTT 9 n^n-t^)/ pj^^^ ■,a,b,C;q) < OO. 

{abq/c, bq, g; g)n (1 - abq) (-ac)" 
In order to achieve this we note that since (or/: q)^ = {aq: q)k{l — aq^~^^)/{l — aq), we have 

c,. 2 , X A 1 - ag'^+i (g~";^)fc(a6g"+i;g)fc(ag;g)fc q^ 

Pn{aq ■,a,b,c-q) = > — ^7 ^ 7 ^ 

1-aq {aq;q)k{cq;q)k {Q;Q)k 

^ 1 \^ {Q^'';Q)k{abq''+^;q)k{aq;q)k q^ x-ipr u \ 
V 7 r-f ^ 7 ^ = {^-aq) Pn{aq;a,b,c;q). 

{aq;q)k{cq;q)k {q;q)k 

Therefore, the series for {4'aq^^'^aq^) is majorized (up to the finite constant (1 — aq)~^) by the 
corresponding series for {ipaq,'^ag)- Thus, ■ijjaq2{x) is an cigcnfunction of I2 and the point ag^ 
belongs to the spectrum of the operator l2- Setting A = aq'^ in (4.7) and acting similarly, one 
obtains that ip^qS (x) is an eigenfunction of I2 and the point ag^ belongs to the spectrum of 
l2- Repeating this procedure, one sees that tpaq^ix), n = 1, 2, • • •, are eigenfunctions of I2 and 
the set og", n = 1, 2, • ■ ■, belongs to the spectrum of I2. Likewise, one concludes that ijjcqn-{x), 
n = 1, 2, • • •, are eigenfunctions of I2 and the set eg", n = 1, 2, • • •, belongs to the spectrum of 
l2- Note that so far we do not know whether the operator I2 has other spectral points or not. 
In order to solve this problem we shall proceed as in subsection 3.2. 

The functions tpaq'^ix) and tpcq^ix), n = 1, 2, ■ ■ ■, are linearly independent elements of the 
Hilbert space TC. Suppose that ag" and eg", n = 1,2, • • •, constitute the whole spectrum of 
the operator I2. Then the set of functions tpaq^ix) and ipcq^ix), n = 1,2,---, is a basis of the 
space H. Introducing the notations := Cag"+i(a^) and := icq"+^{x), = 0, 1, 2, • • •, we 
find from (4.7) that 

JH„ = a-^eg-2"-i(l-ag"+i)(l-6ag"+Vc)H„+i + d„S„ + a-ieg-2"(l-g")(l-ag7e)S„_i, 

JS; = c-iag-2"-i(l - cg"+i)(l - 6g"+i) + < S„ + c-\q-^^{l - g")(l - cg"/a) 
where 

d„ = i[g-2"-ic(l + g) - g-"(a6 + ac + a + c)\, 

d'n = ^te~^"^"^«(l + 9) -g~"(a6 + ac + a + e)]. 

As we see, the matrix of the operator J in the basis S„ = ^aq"+"^{x), = icq"+'^{x)i 
n = 0, 1,2, is not symmetric, although in the initial basis /„, n = 0, 1,2,---, it was 
symmetric. The reason is that the matrix M := {{amn)mn=o (^mn)mn=o) ^i^'^ entries 

amn ■= Pm{aq'^), a'mn ■= Pm{cq'^), m, n = 0, 1, 2, • • • , 

where Pmidq^), d = a,c, are coefficients (4.4) in the expansion ipdq"{x) = J2m (^midq^) fn{x) 
(see above), is not unitary. This matrix M is formed by adding the columns of the matrix 
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{a'^n) columns of the matrix (a^„) from the right, that is, 



M 



V 



ail • 
021 • 


• • aik ■ 

• • a2k ■ 


• 4i • 


• «iz • • • 

• «2Z • • • 


aji ■ 


■ ■ O-jk ■ 







It maps the basis {/n} into the basis {^/'agn+i , ^^^n+i} in the Hilbcrt space Ti. The nonunitarity 
of the matrix M is equivalent to the statement that the basis S„ := ^aqn+i(a;), H„ := ■^cq"+i (2^)) 
n = 0, 1, 2, • • •, is not normalized. In order to normalize it we have to multiply H„ by appro- 
priate numbers c„ and by numbers c^. Let S„ = cViS„, = c^S„, n = 0, 1, 2, • • •, be a 
normalized basis. Then the operator J is symmetric in this basis and has the form 



+c-lic„a-icg-2"(l - aq^/c){l - q^) H„_i, 



; = c';iic'„c-iag-2-i(l - - cq^+^) + S„ 



+c';!ic'„c-iag-2"(l - cg7a)(l - g") S„_i, 
The symmetricity of the matrix of the operator J in the basis {E]„, S^^} means that 



(4.11) 
(4.12) 



c;|ic„g"2"-^(l - a(?"+^)(l - ahq^+'/c) 



c„ic„+ig-2"-2(l - ag"+Vc)(l - 



/-I „-2n-l(i _ ^,gn+l-)(i _ ^ c';^c'„+ig-2"-2(l - C9"+Va)(l - g"+^). 



that is, 

G 



(1 - ag")(l - abq^'/c) 



Cn-l 



Thus, 



Cn — C" 



(l-g")(l-ag"/c) ' <_i 

niabq/c, aq;q)nY^'^ 



, c'^ = C' q 



^ (1 - cg'^)(l - bq'') 
'(l-g")(l-cg"/a) 

t (bq,cq;q)n Y^"^ 
(cq/a, q; q\ 



{aq/c,q;q\ 

where C and C" are some constants. 
Therefore, in the expansions 

1paqn{x) = En{x) = ^ C„ /3^(og'') /^(x) = ^ amnfLi^), 
m m 

Ipcqnix) = En{x) = C'^l3m{cq'')fln{x) = a'^^fl,{x), 

the matrix M := {{a. 

amn = Cn Pmiaq"') = C [q 



(4.13) 
(4.14) 



m 



with entries 



^nJrn,n=0 \"'mn>m,n=0> 

[abq/c, aq; q)n {abq, aq, cq; q)m (1 - abq^'^+^ ] 



iaq/c,q;q)n {abq/c,bq,q;q)m (1 - abq){-ac) 
xq-^("'+^y^Pm{aq''+';a,b,c;q), 



1/2 



(4.15) 



O-'mn = CnPmicq'') = C { q 



1 /9 

{bq, cq; q)n {abq, aq, cq; q)m (1 - abq'^'^+^) \ 



{{cq/a, q; q)n {abq/c, bq, q; q)m (1 - abq){-acy 
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X q-m{m+z)/i .a,h,c;q), (4. 16) 

is unitary, provided that the constants C and C are appropriately chosen. In order to calculate 
these constants, one can use the relations 

m=0 m=0 

for n = 0. Then these sums are multiples of the sums in (4.9) and (4.10), so we find that 

_ {bq, eg; qfj,^ , _ {aq, abq/c; qfj,^ 

[ahq^,cla\q)^ [ahq^ .aj c\q)^ 

The coefficients c„ and in (4.13)-(4.16) are thus real and equal to 

^ ^ ( ibq,cq;q)^ {abq/c, aq; q)n q'' \^^'^ ^ ^ f {aq, abq/c; q)^^ [bq, eg] q)nq'' V'' 
\{abq^,c/a;g)oo {aq/c,q;q)n J ' " \{abq'^,a/c;q)oo icq/a,q;q)n J 

The orthogonality of the matrix M = ((amn)^,n=o i^'mn)m,n=o) means that 



m m m 



^ ], {0'mnO'm'n "I" ^mn^m'n) — ^mm' ■ (4.19) 
n 

Substituting the expressions for amn and into (4.19), one obtains the relation 

{bq, eg; q)^ y. {ag, abg/c; l)n<f p(^^n+l^p^,^^n+l^ 
{abg^,c/a;g)oo {ag/c,g;g)n 



, {aq, abq/c; q)oo ^ {bq, eg; q)ng'^ ^ (,„n+Up r,„n+u 
{abq^,a/c;g)oo ^ {cq/a,q;q)n 



_ {l-abq){bg,abq/c,q;q)m ^ „„^m„m(m+3)/2 . 
- (1 - abq^rn+i^{aq, abq, cq; q^ ^ ^ ^ 



(4.20) 



This identity must give an orthogonality relation for the big g-Jacobi polynomials Pm{y) = 
Pm{y] CL, b, c; q). An only gap, which appears here, is the following. We have assumed that the 
points ag" and eg", ra = 0, 1, 2, • • •, exhaust the whole spectrum of the operator I2. As in the 
case of the operator /i in subsection 3.2, if the operator I2 had other spectral points Xk, then on 
the left-hand side of (4.20) would appear other summands fixf,Pm{xk] a, b, c; q)P^i{xk; a, b, c; q), 
which correspond to these additional points. Let us show that these additional summands do 
not appear. We set m = m' = in the relation (4.20) with the additional summands. This 
results in the equality 



{bq,cq;q)oo ^ {aq,abq/c;q)nq'^ 

2 r./r,-„\ 



,c/a;q)oo {aq/c,q;q) 

{aq,abq/c;q)oo {bg,cq;q)nq 



{abq'^,a/c;q)oo {cq/a,q;q)n ^ 
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In order to show that ^^i. ii^k = 0) ^Sike into account the relation 

{Aq/C,Bq/C-q), 



{q/C,ABq/C-q) 



°°2MAB;C; q,q) 



+ (r/^'^fnlr ^ MAg/C, Bq/C- q'/C- q,q) = l 

{C/q,ABq/C;q)oo 

(see formula (2.10.13) in [GR]). Putting here A = aq, B = abq/c and C = aq/c, we obtain 
relation (4.21) without the summand fixk- Therefore, in (4.21) the sum fi^^. does really 
vanish and formula (4.20) gives an orthogonality relation for big g-Jacobi polynomials. 

By using the operators I2 and J, we thus derived the orthogonality relation for big g-Jacobi 
polynomials. 

The orthogonality relation (4.20) enables one to formulate the following statement: The 
spectrum of the operator I2 coincides with the set of points aq^'^^ and cq^'^^ , n = 0, 1, 2, • • •. 
The spectrum is simple and has one accumulation point at 0. 

4.3. Dual big g-Jacobi polynomials 

Now we consider the relations (4.18). They give the orthogonality relation for the set of 



a' 

coincide with the functions 



matrix elements Omn ^^<i o!mni viewed as functions of m. Up to multiplicative factors, they 



F„(x;a,6,c;g) := 302(x,a6g/x,ag"+^; ag,cg; g,g), n = 0, 1, 2, • • • , (4.22) 

F^(x;a, 6,c;g) := 3(/)2(x,a6g/x,cg"^+^ ag,cg; q,q) = Fnix;c,ab/c,a), n = 0, 1,2, (4.23) 
considered on the corresponding sets of points. Namely, we have 

a =a (a b c) - C ( ("^g/'^^ "g^ g)" ("^g^ "g^ ^9? l)m (1 - Q^g^"""^^) \ 
mn - mn , , ^ {aq/c,q;q)n {abq/c,bq,q;q)m{^ - abq){-ac)"' J 

amn-amn{a,0,C) U ^q (^^/^^ ^. (^^^/^^ ^. (1 _ (_„^)m ^ 

xq-"'^^+'y^F;,{q-^;a,b,c;q) = amn{c,ab/c,a), (4.25) 

where C and C are given by formulas (4.17). The relations (4.18) lead to the following 
orthogonality relations for the functions (4.22) and (4.23): 

(6g,cg;g)oo ^ n 77 / -m u \ i? 1 -m u \ {m/c,q;q)n <■ . , „«x 

— ^ — -, ^ ^ p{m)Fn{q ■,a,b,c;q)Fn,{q "'■a,b,c;q) = — r — -Sn„', (4.26) 

{abq^,c/a;q)oo {aq,abq/c;q)nq" 

{aq, abq/c; g)oo ^ , , , n (cg/a, g; g)n . 

p^,^^^^ X; p(m)F„(g ;a,6,c;g)F„,(g ; a, 6, c; g) = ^^^^^^^ (4.27) 

00 

p{m)Fn{q~"';a,b,c;q)F;^,{q-^;a,b,c;q)=0, (4.28) 

m=0 

where 

/ ) (1 - a&g^"'+^)(ag, abq, cq; q)m ^-m(m+3)/2_ 
(1 - abq){bq, abq/c, g; q)m (-ac)'" 
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There is another form for the functions Fn{q "^;a,b,c;q) and Fl^{q ^\a,h,c\q). Indeed, 
one can use the relation (III. 12) of Appendix III in [GR] to obtain that 



Fn{q "'■,a,b,c;q) = 



{cq '"MgV ^^^^^^i^^ 
(eg; q)m 



q~"',abq"'+'^,q-'' 
aq, abq/c 



q,aq^+'/c 



{abq/c; q)r 



-"',abq"'+^,q-'' 



q,aq^+'/c 



and 



{aq;q)m 



F'^{q-^;aAc;q) 



aq, abq/c 

q-"',abq"'+^,q~'' 
bq, cq 



q,cq^+'/a]. 



The basic hypergeometric functions 3(^)2 in these formulas are polynomials in /i(m) := g ™ + 



So if we introduce the notation 
Dn{lJ,{m);a,b,c\q) := 3(j)2 



q-"\abq"'+^,q-'' 
aq, abq/c 



q,aq^+^/c], 



(4.29) 



then 



F„(g-;a,6,c;g) = (-c)-g-(-+^)/^i^.(/x(m); a, 6, c|g), 

(,cg, qjm 

Fl,{q-^;aAc;q) = p^{-arq^^^+^y''DMm)-Aa,ab/c\q). 
[aq, qjm 

Formula (4.26) directly leads to the orthogonality relation for the polynomials Dn{ji{m)) = 
Dn{^x{m);a,b,c\q): 



E 

m=0 



(1 — abq^'^^^){aq, abq, abq/c; q)r 



(1 - abq){bq,cq,q;q)m 

{abq^, c/ a; q) 00 {aq/c,q;q) 



i-c/ar L>4/x(m)) Dn,{,,{m)) 



ibq,cq;qU (aq, abq/ c; qW^'^' ' ^''''^ 

From (4.27) one obtains the orthogonality relation for the polynomials Dn{iJ,{m);b,a,ab/c\q) 
(which follows also from the relation (4.30) by interchanging a and b and replacing c by ab/c). 

We call the polynomials Dn{^{m);a,b,c\q) dual big q-Jacobi polynomials. It is natural to 
ask whether they can be identified with some known and thoroughly studied set of polynomials. 
The answer is: they can be obtained from the g-Racah polynomials Rn{lJ^{x);a,b,c,d\q) of 
Askey and Wilson [AW] by setting a = q~^~^ and sending N ^ 00, that is, 



Dnin{x);a,b,c\q) = liia Rn{n{x);q ^ ^,a/c,a,b\q). 

iV— >oo 



(4.31) 



Observe that the orthogonality relation (4.30) can be also derived from formula (4.16) in 
[Ros]. But the derivation of this formula (4.16) is rather complicated. 

The dual polynomials (4.29) and the big g-Jacobi polynomials (4.3) are interrelated in the 
following way: 



DnifJ'im);a,b,c\q) 



(eg; g)r 



{abq/c; q\ 
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It is worth noting here that in the hmit as c — > the dual big g-Jacobi polynomials 
Dn{fJ,{x);a,b,c\q) coincide with the dual little q-Jacobi polynomials dn{fJ'{x);b,a\q), defined 
in section 3. The dual little g-Jacobi polynomials dn{n{x);a,b\q) reduce, in turn, to the Al- 
Salam-Carlitz II polynomials Vn'^\s;q) on the g-linear lattice s = q~^ (see [KSw], p. 114) in 
the case when the parameter b vanishes, that is, 

cZ„(/x(ar);a,0|g) =2^0(9-", g"";-; g,g7«) = (-a)-^"^"-^)/^ g). (4.32) 

This means that we have now a complete chain of reductions 

Rn{n{x);a,b,c,d\q) — > Dn{l^{xy,c,d,c/b\q) — y dn{fi{xy,d,c\q) —^Vji'^\q~'-';q) 

a—*oo 0— >oo c=U 

from the four-parameter family of g-Racah polynomials, which occupy the upper level in the 

Askey-scheme of basic hyper geometric polynomials (sec [KSw], p. 62), down to the one- 
parameter set of Al-Salam-Carlitz II polynomials from the second level in the same scheme. 
So, the dual big and dual little g-Jacobi polynomials L'„(/z(x); a, 6, c|g) and (i„(/z(x); a, 6|g) 
should occupy the fourth and third level in the Askey-scheme, respectively. 

The recurrence relations for the polynomials I?„(/x(m) = Z)„(/x(m); a, 6, c|g) are obtained 
from the g-difference equation (4.6). It has the form 

^q-m _ _ abq^+')DMm)) = A£'n+i(/x(m)) - {A^ + C„)L>„(/x(m)) + C„Z?„_i(/x(m)), 
where 

An = g-'"-Hl - a^"^') [{c/a) - bq^+'] , C„ = ^-^"(l - q^) [(c/a) - g"] . 

The relation (4.28) leads to the equality (another proof of this relation is given in Ap- 
pendix) 

/ -.m — cbq'^"^^^)(abq:q)m "<-(m-i) „ , , , , , ^ ^ , , ^ , , , , n 
(~^) n _ hV \ DMmy,a, b, c\q) Dn'{fi{my, b, a, ab/c\q) = 0. 

(4.33) 

Note that from the expression (4.29) for the dual big g-Jacobi polynomials it follows that 
they possess the symmetry property 

-D„(/i(m); a, b, c\q) = Dn{n{myab/c, c, b\q). (4.34) 

The set of functions (4.22) and (4.23) form an orthogonal basis in the Hilbert space P of 
functions, defined on the set of points m = 0, 1, 2, • • •, with the scalar product 

00 

(/i, /2) = Yl /2M' 

m=0 

where p(m) is the same as in formulas (4.26)-(4.28). Consequent from this fact, one can 
deduce (in the same way as in the case of dual little g-Jacobi polynomials) that the dual big 
q-Jacobi polynomials DnifJ-im); a, b, c\q) correspond to indeterminate moment problem and the 
orthogonality measure for them, given by formula (4-30), is not extremal. 
It is difficult to find extremal measures for these polynomials. 
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4.4. Generating functions for dual big ^-Jacobi polynomials 

Generating functions are known to be of great importance in the theory of orthogonal 
polynomials (see, for example, [AAR] and [SM]). For the sake of completeness, we briefly 
discuss in this section some instances of linear generating functions for the dual g-Jacobi 
polynomials Dn{fi{x); a, b, c\q) and d„(/Lt(x); a, b\q). To start with, let us consider a generating- 
function formula 

Y,(^eDM.)^aM,) = ,<h i'"','^" 1, mtic) , (4.35) 

^ {(l\(l)n {t;q)oo V (^bq/c, aqt ) 

where \t\ < 1 and, as before, //(x) = q~^ + abq^~^^. To verify (4.35), insert the explicit form 
(4.29) of the dual big g-Jacobi polynomials 



r> f f \ u \ \ '',abq''+\q ";gfc /ag"+^ 

r»n(/x(x); a, b, c\q) = > — — 

^ {aq,abq c,q;q)k \ c 



k=o 



k 



into the left side of (4.35) and interchange the order of summation. The subsequent use of 
the relations 

{a;q)m+k = {a;q)m{aq"';q)k = {a; q)k{aq''; q)m, 
{q-"'-'';q)k = {-ifq-'^^-^^^+'^^'^q^^^^qh 

(sec [GR], Appendix I) simplifies the inner sum and enables one to evaluate it by the q- 
binomial formula (3.16). This gives the quotient of two infinite products in front of 2^2 on 
the right side of (4.35), times {aqt] q)^^- The remaining sum over k yields 2(^2 series itself. 

As a consistency check, one may also obtain (4.35) directly from the generating function 
for the g-Racah polynomials Rn{iJ,{x);a,f3,^,6\q) (see formula (3.2.13) in [KSw]) by setting 
a = q~'^~^ and sending 00. This results in the relation 

W1)n (i;g)oo V o-bq/c ) 

The left side of (4.36) depends on the variable x by dint of the combination /x(x) = q"^ + 
ab(f^^. Off hand, it is not evident that the right side of (4.36) is also a function of the 
lattice ii{x). Nevertheless, this is the case. Moreover, the right sides of (4.35) and (4.36) are 
equivalent: this fact is known in the theory of special functions as Jackson's transformation 

2<?!)i(a,6; c; q,z) = 24>2{a,c/b; c,az; q,bz) 

Qjoo 

(see, for example, [GR]). 

The symmetry property (4.34) of the dual big g-Jacobi polynomials Dn{n{x);a,b,c\q), 
combined with (4.35), generates another relation 

{abq/c;q)n,n^ , , . , , x {abqt/c;q)^ f q'"", abq''+^ \ 

— r eDn{fJ,{x);a, b, c\q) = — — 2(p2 [ , , , q, aqt c 

{Q;Q)n {t;q)oo V abqt c J 



{abtq^^'/c-qU . ( <r\ h'^-^ ^^ ^^^ 
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Similarly, a generating function for the dual little g-Jacobi polynomials has the form 

One can verify (4.38) directly by inserting the explicit form (3.20) of d„(//(x); a, 6|g) into the 
left side of (4.38) and repeating the same steps as in the case of deriving (4.35). This will lead 
to the expression 

( nhni-' n\-- 

■20i(g"^, abq^'^^; abqt; q,t) 



{at; q)c 

and it remains only to employ Heine's summation formula (1.5.1) from [GR]. After a simple re- 
scaling of the parameters the generating function (4.38) coincides with that, obtained earlier 
in [BI]. 

The simplest way of obtaining (4.38) is to send c — in both sides of (4.36): the 2(l>i series 
on the right side of (4.36) reduces to i(l>o{q~^; — ; q, t/a), which is evaluated by the g-binomial 
formula (3.16). 

Finally, when the parameter b vanishes, (4.38) reduces to the known generating function 



^ {q;q)n " '^^ {at,t;q) 



oo 
oo 



n=0 

for the Al-Salam-Carlitz II polynomials (see (3.25.11) in [KSw]). 

5. DISCRETE g-ULTRASPHERICAL POLYNOMIALS AND THEIR DUALS 
5.1. Discrete g-ultraspherical polynomials 

For the big g-Jacobi polynomials Pn{x;a,b,c;q) the following limit relation holds: 

limP (x-a°' -a^-a) = — — 

limt^^[x,q ,q , q ,q) p(a,0) ^-^y 

where 7 is real. Therefore, limqji Pn{x; q", q°', —q^'] q) is a multiple of the Gcgenbauer (ultra- 
spherical) polynomial '^^'^\x). For this reason, we introduce the notation 

ClT!''\x;q) := Pn{x;a,a,-a;q) = 3(?!>2(g"", a^g'^+S x; aq,-aq; q,q). (5.1) 

It is obvious from (5.1) that Cn\x;q) is a rational function in the parameter a. 

From the recurrence relation for the big g-Jacobi polynomials (see subsection 4.1) one 
readily verifies that the polynomials (5.1) satisfy the following three-term recurrence relation: 

xC(«)(x;g) = ^„(a)ci^i(x;g) + C„(a)ci"i(x;g), (5.2) 

where An{a) = (1 - ag"+i)/(l - ag2"+i), C7„(o) = 1 - ^„(a), and C^''\x;q) = 1. 

An orthogonality relation for Cn\x]q), which follows from that for the big g-Jacobi 
polynomials and is considered in the next section, holds for positive values of a. We shall see 
that the polynomials cl^\x; q) are orthogonal also for imaginary values of a and x. In order 
to dispense with imaginary numbers in this case, we introduce the following notation: 

Clf\x;q) := (-i)"C(-«')(ix; g) = (-i)"P„(irr; ia, ia, -ia; g), (5.3) 
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where x is real and < a < cx). The polynomials cii'\x; q) satisfy the recurrence relation 

xCi-\x; q) = An{a) cl^l^{x- q) + (7„(a) cl^l^{x- q), (5.4) 

where i„(a) = An{-a) = (1 + ag"+i)/(l + ag^^+i), Cn{a) = in(a) - 1, and C^''\x;q) = 1. 
Note that An{a) > 1 and, hence, coefficients in the recurrence relation (5.4) for cii'\x;q) sat- 
isfy the conditions An{a)Cn+i{a) > of Favard's characterization theorem for n = 0, 1, 2, • • • 
(see, for example, [GR]). This means that these polynomials are orthogonal with respect to a 
positive measure with infinitely many points of support. An explicit form of this measure is 
derived in the next section. 
So, we have 



q,q]. (5.5) 



Ci^ix;q) ^ i-^-C^^i^x;q) ^ (-^^^.^.^^^ 

(Here and everywhere below under ^/a, a > 0, we understand a positive value of the root.) 
From the recurrence relation (5.4) it follows that the polynomials (5.5) are real for x G M and 
< a < DO. From (5.5) it is also obvious that they are rational functions in the parameter 
a. Observe that the situation when along with orthogonal polynomials Pn{x), depending on 
some parameters, the set of polynomials (— i)"p„(ix) is also orthogonal (but for other values 
of parameters) is known; see, for example, [Rom], [Ask], and [CS]. 

We show below that the polynomials Cn\x;q) and Cn\x;q), interrelated by (5.5), are 
orthogonal with respect to discrete measures. For this reason, they may be regarded [AK4] 
as a discrete version of g-ultraspherical polynomials of Rogers (see, for example, [AI2]). 

Proposition 5.1. The following expressions for the discrete q-ultraspherical polynomials 
(5.1) hold: 

= ^^^4^i-^)'^'^'^'W/aq'-,q-\a\q% (5.6) 



Citi(^;q) = ^f''f\f {-l)'q'^'+'^xp,{x'/aq'; q, a\q% (5.7) 
[aq^'-iQ )k 

where pk{y;a,b\q) are the little q-Jacobi polynomials (3.3). 

Proof. To start with (5.6), apply Singh's quadratic transformation for a terminating 3^2 
series 

' ' ' ^ ' g , g , (5.8) 



3^2 



abq^/"^, —abq^/ 



392 



a'^b'^q, 



which is valid when both sides in (5.8) terminate (see [GR], formula (3.10.13)), to the expres- 
sion in (5.1) for g-ultraspherical polynomials C^{x;q). This results in the following: 

C^?(x;,) =3</>2 aq'^^\ x'- aq\ 0; q^) . 

Now apply to this basic hypergeometric series ■1,^2 the transformation formula 

b 

c 

(see formula (III.7) from Appendix III in [GR]) in order to get 



29^1 



{c/b-q)n , 6, bzq-^/c 



q,q\ (5.9) 
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Comparing this formula with the expression for the httle g-Jacobi polynomials (3.3) one arrives 
at (5.6). 

One can now prove (5.7) by induction with the aid of (5.9) and the recurrence relation 
(5.2). Indeed, since Cg"'*(,T;q) = 1 and Ao(a) = 1, one obtains from (5.2) that c["'\x;q) = x. 
As the next step use the fact that c'^\x\q) = ^(t)2{q~^ ,aq^ ,x^; aq^,0; q^,q^) to evaluate 
from (5.2) explicitly that 

C'^\x\q) = x^(l)2{q~'^, aq^, x^; aq^, 0; 
So, let us suppose that 

Ct^_,ix;q) = X3<^2 {q-'^'-'\ aq''^\ x'; aq\ 0; q\q') (5.10) 

for = 1,2,3,---, and evaluate a sum A2;3(a)xC2^^ (a:;; g) + (l — ^2fc (a))C2^^_j^(x; g). As follows 

from the recurrence relation (5.2), this sum should be equal to c'^^j^^{x; q) . This is the case 
because it is equal to 



q ag^^+^, x^ 
aq^, 



q\qA+(\-A-})^^2 



ag^, 



2 2 

q ,q 



= a; 302 



q aq^^^^, x^ 



aq^, 

if one takes into account readily verified identities 



2 2 

q ,q 



(5.11) 



^2fc(5 ^'';9^)m -r VJ- - --^2fc 

1 _ ag2(fe+m)+l 



+ {l-A,')iq-'^'-'^;q')m = K^ 



3{k+m)+l 



-{q-"';q') 



1 _ a^2fc+i 

2,^, {aq''''-';q')m={aq"'^';q')m, 



1 — aq 

for m = 0, 1, 2, • • • , k. The right side of (5.11) does coincide with C'2^^_j_j(aj; g), defined by the 
same expression (5.10) with k — > k+1. Thus, it remains only to apply the same transformation 
formula (5.9) in order to arrive at (5.7). Proposition is proved. 

Remark. Observe that in the process of proving formula (5.7), we established a quadratic 
transformation 



3^2 



q "^^ ^,aq^^^'^,x 
y/aq, -y/aq 



X3<p2 



q "^^ , aq^^^^ , x^ 
aq"^, 



2 2 

q ,q 



(5.12) 



for the terminating basic hypergeometric polynomials 3^2 with k = 0, 1, 2, • • •. The left side 
in (5.12) defines the polynomials C^^^iix; q) by (5.1), whereas the right side follows from the 
expression (5.11) for the same polynomials. The formula (5.12) represents an extension of 
Singh's quadratic transformation (5.8) to the case when = g-^fe-i ^^^^^ therefore, the left 
side in (5.8) terminates, but the right side does not. 

It follows from (5.5)-(5.7) that 



cS{x;q) = " i-l^q'^'^^UxVaq'; q-\ -a\q% 



(5.13) 
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Citii^;<l) = f''i\i {-l)V^'^''^xp,{x'/aq';q, -a\q'). (5.14) 

In particular, it is clear from these formulas that the pol3momials Cn\x;q) are real-valued 
for X G M and a > 0. 

5.2. Orthogonality relations for discrete g-ultraspherical polynomials 

Since the polynomials C^\x; q) are a particular case of the big (7-Jacobi polynomials, an 
orthogonality relation for them follows from (4.20). Setting a = b = — c, a > 0, into (4.20) 
and considering the case when m = 2k and m' = 2k', one verifies that two sums on the left of 
(4.20) coincide (since ab/c = —a = c) and we obtain the following orthogonality relation for 
Ciiix-q): 

00 / 2 2\ s 

E (J'J:' CiiiV^q^^^;q)C^i,iV^q^^^;q) 

(ag3; g^)^ {l-aq)a'^'' {q; q)2k g'^^^fc+s) 



iq;q^U l-ag^^^+i {aq;q) 



2k 



^kk' 



(5.15) 



where ^/a, a > 0, denotes a positive value of the root. Thus, the family of polynomials 
C^j}{x;q), k = 0, 1,2, with < a < q~^, is orthogonal on the set of points ^/aq'^~^^, 
s = 0,l,2,---. 

As we know, the polynomials C^{x;q) are functions in x^, that is, Cfj{^q'+^;q) IS m 

fact a function in ag^*+^. The set of functions Cg^^ (x; g). A; = 0, 1, 2, • • •, constitutes a complete 
basis in the Hilbert space of functions f{x^) with the scalar product 

s=0 W ,9 j5 

This result can be obtained from the orthogonality relation for the little qf-Jacobi polynomials, 
if one takes into account formula (5.6). 

Putting a = b = — c, a > 0, into (4.20) and considering the case when m = 2A; + 1 and 
m! = 2k' + 1, one verifies that two sums on the left of (4.20) again coincide and we obtain the 
following orthogonality relation for C2^^_j_^(x; g): 



/ 2 2\ s 



s=0 ^« 



_ (Qg^; g^)oo (1 - aq) (g; g)2fc+l „(fc+2)(2fc+l) , 

~ iq;q')oo {l-aq'^+^) {aq;q)2k+,^ ^ ^ 

The polynomials C'2^^^^(x; q), A; = 0, 1, 2, • • with < a < q""^ , are thus orthogonal on the set 
of points y/aq^'^^ , s = 0, 1, 2, • • •. 

The polynomials x~^C^j^'^{x;q) are functions in x^, that is, x~^ C^j}j^'^{y/a q'^'^^ ; q) are in 

fact functions in aq^^^"^ . The collection of functions C^^tj^ii^'^ q) ^ k = 0, 1,2, • • •, constitute a 
complete basis in the Hilbert space of functions of the form F{x) = x/(x^) with the scalar 
product 

00 / 2 2\ s 

(i^„i^2) 72%^ F,{V^q-^+')F2{V^q^^^). 

s=0 W >g 
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Again, this result can be obtained from the orthogonahty relation for the little g-Jacobi 
polynomials if one takes into account formula (5.7). 

We have shown that the polynomials C^'^^x; g), A; = 0, 1, 2, • • •, as well as the polynomi- 
als C2^^^(a;;g), k = 0, 1,2, • • •, are orthogonal on the set of points ^/aq^'^^, s = 0, 1,2, • • •. 

However, the polynomials c'!^{x;q), k = 0, 1,2, are not orthogonal to the polynomials 

C^^j^^{x; q), = 0, 1, 2, • • •, on this set of points. In order to obtain an orthogonality for the 

whole collection of the polynomials Cn"'' {x',q), n = Q,\,2, - • • , one has to consider them on the 
set of points ±^/aq^^^, s = 0, 1,2, • • •. Since the polynomials from the first set are even and 
the polynomials from the second set are odd, for each A;, A;' G {0, 1, 2, ■ ■ ■} the infinite sum 

coincides with the following one 

oo ^ 2 2\ s 

T — ]M^Q_jsQ_ Ma) ( /t: „s+l . „\ r<{a) ( f7:„s+l.„\ 



s=0 



Therefore, I1 — I2 = 0. This gives the orthogonality of polynomials from the set {x; q), k = 
0, 1, 2, • • •, with the polynomials from the set C'2fc+i('^' 'S')' ^ ~ 0) !> 2, • • •. The orthogonality 
relation for the whole set of polynomials Cn\x; g), = 0, 1, 2, • • •, can be written in the form 

00 / 2 2\ s 

s=Oe=±l W >9 

_ iaq';q^)oo (1 - ag) a" (g; g)n g"("+^)/% 

- {q;q^U (1 - ag^n+i) (aq;q)n ^'-'"^ 

We thus see that the polynomials Cn\x; q), n = 0,1,2, ■■■ , with < a < q~^ are orthogonal 
on the set of points ±^/aq'^~^^, s = 0, 1, 2, • • •. 

An orthogonality relation for the polynomials Cn\x; g), n = 0, 1, 2, • • •, is derived by using 
the relations (5.13), (5.14), and the orthogonality relation for the little g-Jacobi polynomials. 

Writing down the orthogonality relation (3.15) for the polynomials pfc(.T^/ag^; g~^, —ajg^) 
and using the relation (5.13), one finds an orthogonality relation for the set of polynomials 
C^^{x; q), A: = 0, 1, 2, • • •, with a > 0. It has the form 

00 / 2 2\ s 

E ^^^^ Ct^^aq^^^q) C^f, (V^ q^^^ ; q) 



s=0 ^« 



{-aq : q )^ (1 + uq) a {q: q)2k M2k+3) , ...on 

= ~? 2^ TT"^ ik+T\7 \~1 (5.18) 

(g; g^)oo (1 + ag^'^+i) {-aq; q)2k 

Consequently, the family of polynomials c'^^{x; q), = 0, 1, 2, • • is orthogonal on the set of 
points Y^g*^^; s = 0, 1, 2, • • • . 

As in the case of polynomials c'^f}{x; q), k = 0, 1, 2, • • •, the set c'i^i^ {x; q), A; = 0, 1, 2, • • •, 
is complete in the Hilbert space of functions /(x^) with the corresponding scalar product. 
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Similarly, using formula (5.14) and the orthogonality relation for the little g-Jacobi poly- 
nomials Pk{x^ / o-ff 'i W6 fi^d the orthogonality relation 



s=o Is 



_ (-ag^^ q^)^ (1 + aq) a?^+^ {q; q)2k+i jk+2)(2k+i) , 1Q^ 
" {q;Q')oo {l + aq^^^^) {-aq-q),k+,^ '''' ^'-'"^ 

for the set of polynomials C2'k+i(^'^ q), k = 0,1,2,- ■ ■. We see from this relation that for a > 

the polynomials (52^+1 k = 0,1,2, ■■■ , are orthogonal on the same set of points ^/aq'^~^^, 

s = 0,l,2,---. 

Thus, the polynomials C'2^'*(x; (7), k = 0,1,2,- ••, as well as the polynomials C2^^^{x;q), 
k = 0,1,2,- ■■, are orthogonal on the set of points ^/aq'^^^, s = 0, 1,2,---. However, the 

polynomials C'2fe^(^' k = 0,1,2,- ■ ■, are not orthogonal to the polynomials C'2fc+i(^' l)i ^ = 
0, 1, 2, • • • , on this set of points. As in the previous case, in order to prove that the polynomials 

^2ki^'i Q)^ k = 0,1,2,- ■ -, are orthogonal to the polynomials C'2fc+i(^' q): k = 0,1,2,- - -, one 
has to consider them on the set of points ±^/aq^^^, s = 0, 1,2, • • •. Since the polynomials 
from the first set are even and the polynomials from the second set are odd, then the infinite 
sum 

00 / 2 2\ s 

_s-^ [-aq ]q )sQ Ma)/ z;: r^(") ( /7;nS+^-n\ 



coincides with the sum 



00 



i2 = -2^ — / 2.^2) ^2kK-VaQ ,gj<-2fc'+i(-vag 

Consequently, Ii — I2 = 0. This gives the mutual orthogonality of the polynomials C^{x; q), 
A; = 0, 1, 2, • • •, to the polynomials C'2fe+i(^' ^ = 0,1,2,- - -. Thus, the orthogonality relation 
for the whole set of polynomials Cn\x; q), n = 0,1,2,- - -, can be written in the form 

E E i2'!J:' C^\eV-aq^^\q)C^}[eV-aq^^\q) 

_ (-ag^;g^)oo (l + ag)a" {q;q)n „n(n+3)/2 . 

" (9; 5^)00 (l + ag2n+i)(_a5;^)„^ ^^"^^^ 

Note that the family of polynomials Cn\x; q), n = 0,1,2, - - -, corresponds to the determinate 
moment problem, since the set of orthogonality is bounded. Thus, the orthogonality measure 
in (5.20) is unique. 

In fact, formula (5.20) extends the orthogonality relation for the big g-Jacobi polynomials 
Pn{x; a, a, —a; q) to a new domain of values of the parameter a. 

5.3. Dual discrete g-ultraspherical polynomials 

The polynomials (4.29) are dual to the big Q'-Jacobi polynomials (4.3). Let us set a = 
b = — c in the polynomials (4.29), as we made before in the polynomials (4.3). This gives the 
polynomials 



s+l. 



\lJ'{x;a^)\q) ■-= Dn{lJ.{x;a'^);a,a,-a\q) := 3(j)2{ 



g,-g"+M, (5.21) 
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where fj,{x; a^) = q ^ + a?q'^~^^. They satisfy the three-term recurrence relation 
{q--+aq-+^)D(-\,,{x-a)\q) = -q-^--\l - ag2"+2)L>(«|,(/,(x; a)|g) 



+q 



-2n-l 



{l+q)D(:\^,{x■,a)\q) - q-'^^l q'-)Dl:lM^; <^M, 



which follows from the recurrence relation for the polynomials Dn{lJ'{x;ab);a,b,c\q) from 
section 4.3. 

For the polynomials Dn ' {fi{x; a^)\q) with imaginary a we introduce the notation 



D^^\lj,(x; -a'^)\q) := Dn{fi{x] -a^);ia,ia, ~ia\q) 
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q ^, —a?q^~^^,q 



laq, —laq 



q,-Q 



(5.22) 



The polynomials D^\ii{x; — a^)|g) satisfy the recurrence relation 

{q-- - ag-+i)^(<^)(//(x; -a)\q) = -q-^^-\l + aq^^+^)b^:l,{p.{x- -a)\q) 

+g-2-i(l + g)^(«)(/x(x; -a)\q) - g-2"(l - g2n)^(«)^(^(^. _«)|^). 

It is obvious from this relation that the polynomials Dn\lJ.{x\ —a)\q) are real for x G M and 
a > 0. For a > these polynomials satisfy the conditions of Favard's theorem and, therefore, 
they are orthogonal. 

Proposition 5.2. The following expressions for the dual discrete q-ultraspherical polyno- 
mials (5.21) hold: 



L>i«)(/i(2fe; a)\q) = d„(/x(fc; q-\)- q-\ a\q^) 
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q , uq , q 

aq^ 



L»(")(/x(2A; + 1; a)\q) = g"d„(/x(A;; ga); g, a\q^) 
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q , , y 



(5.23) 



(5.24) 



where k are nonnegative integers and dni^lJ-ix; be); b, c\q) are the dual little q-Jacobi polynomials 
(3.20). 

Proof. Applying to the right side of (5.21) the formula (III. 13) from Appendix III in [GR] 
and then Singh's quadratic relation (5.8) for terminating 3^2 series, after some transformations 
one obtains 



Df\lJi{2k;a^)\q)=a-^\-^^^'^+^\4>2 
Now apply the relation (0.6.26) from [KSw] in order to get 

Df\y^{2k;a')\q) 



2k „2„2fc+l „2„2ji+2 



a^g^, 



2 2 
Q ,Q 



ia^q^;q^)k 



-291 
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Using formula (III. 8) of Appendix III from [GR], one arrives at the expression for the polyno- 

( 2 ^ 

mials Dn (/x(2fc; a^)|(;) in terms of the basic hyper geometric function from (5.23), coinciding 
with (i„(^(A;;g~^a^);g~-^,a^|g^). 

The formula (5.24) is proved in the same way by using the relation (5.12). Proposition is 
proved. 

For the polynomials D^n^ {^^{m\—a)\(i) with nonnegative integers m, we have the expres- 
sions 

DtHK^k; -a)\q) = dMk; -q-^a);q-\ -a\q^) 



q -aq'^'^^^, q 
* —aq'^ 



(5.25) 



Di^\l^{2k+1; -a)\q) = g"d„(/x(A;; -ga); g, -a|g2) 



q 391 



-aq^ 



It is plain from the explicit formulas that the polynomials Dli'\iJ,{m)\q) and Dn\i^{m)\q) are 
rational functions of a. 

5.4. Orthogonality relations for dual discrete g-ultraspherical polynomials 

An example of the orthogonality relation for D^^ \fi{x; a'^)\q) = Dn{fJ-{x;a'^);a,a,—a\q), 
< a < q~^, has been discussed in section 4. However, these polynomials correspond to 
the indeterminate moment problem and, therefore, this orthogonality relation is not unique. 
Let us find other orthogonality relations. In order to derive them we take into account the 
relations (5.23) and (5.24), and the orthogonality relation (3.21) for the dual little g-Jacobi 
polynomials. By means of formula (5.23), we arrive at the following orthogonality relation for 
< a < g"^: 



where /n(2/c) = ij,{2k;a). The relation (5.24) leads to the orthogonality, which can be written 
in the form 

£ g'=(^^+^)i^(»)(M2^ + l)k) Di:Hf^{2k + l)\q) 

^ [i- - aq)[q,q)2k+i 

(ag3;g2)^ (g2;g2)ng" 



{q;q^)oc {aq^;q^)n 
where iJ,{2k + 1) = m(2A; + 1; a) and < a < g"^. 

Thus, we have obtained two orthogonality relations for the polynomials D^\fjL{x; a)\q) , 
< a < q~^, one on the lattice /u(2A;; a) = q""^^ + aq^^'^^, k = 0, 1, 2, • • •, and another on 
the lattice fi(2k + 1; a) = q^"^^^^ + aq"^^^^, k = 0, 1, 2, • • •. The corresponding orthogonality 
measures are extremal since they are extremal for the dual little g-Jacobi polynomials from 
formulas (5.23) and (5.24) (see section 3). 

The polynomials dIi\ijl{x; a)\q) also correspond to the indeterminate moment problem 
and, therefore, they have infinitely many positive orthogonality measures. Some of their 
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orthogonality relations can be derived in the same manner as for the polynomials Dn\iJ,{x)\q) 
by using the connection (5.25) and (5.26) of these polynomials with the dual little g-Jacobi 
polynomials (3.20). The relation (5.25) leads to the orthogonality relation 



(-ag3;g2)^ {q'^-q'^)^q- 



5„ 



where ii{2k) = jj,{2k;—a), and the relation (5.26) gives rise to the orthogonality relation, 
which can be written in the form 



k=o (1 + «9)(5;^)2fc+i 

{-aq^;q'^)^ (g^;g^)ng"" . 

0, 



nn' ) 

n 



where ^{2k + 1) = ii[2k + 1; —a). In both cases, a is any positive number. 

Thus, in the case of the polynomials D'^\ii{x\—a)\q) we also have two orthogonality 
relations. The corresponding orthogonality measures are extremal since they are extremal for 
the dual little g-Jacobi polynomials from formulas (5.25) and (5.26). 

Note that the extremal measures for the polynomials D'>a\^M{x)\q) and D^n\p'{x)\q), dis- 
cussed in this section, can be used for constructing self-adjoint extensions of the closed sym- 
metric operators, connected with the three-term recurrence relations for these polynomials 
and representable in an appropriate basis by a Jacobi matrix (details of such construction are 
given in [Ber], Chapter VII). These operators arc representation operators for discrete series 
representations of the quantum algebra Uq{sni^i) (sec, for example, [KS] for a description of 
this algebra). Moreover, the parameter a for these polynomials is connected with the number 
Z, which characterizes the corresponding representation of the discrete series. 

5.5. Other orthogonality relations 

The polynomials D^\ii{x; a)\q) and the polynomials D^\ii{x;—a)\q) correspond to the 

indeterminate moment problems. For this reason, there exist infinitely many orthogonality re- 
lations for them. Let us derive some set of these relations for the polynomials Dn\^{x; —a)\q), 
by using orthogonality relations for the polynomials (5.18) in [BI]. These polynomials are (up 
to a factor) of the form 

n4(e« - e-^)/2-MMq) = 3^i ( ^'^^*''_~2ji~if '^~1?-9%A2) (5.27) 
and orthogonality relations, parameterized by a number d, q < d < 1, are given by the formula 

(-iig-"/d, -t^q-'^/d, t2q''d- q)^ (^4n^n(2n-l) ^ ^2^2n) 



E 



{-tit2/q; q)oo {-(P;q)ooi-q/(P; q)ooiq; q)c 

{q;q)r{ti/t2) 



xur {{d-^q-- - dq^)l2-MM) Us ((d-^g- - dq^)/2-MM) = ^ /I / \ M (5-28) 

{-q^/tit2;q)rq'^ 
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The orthogonality measure here is positive for ti,t2 and tit2 > 0. It is not known whether 
these measures are extremal or not. 

In order to use this orthogonality relation for the polynomials -D^"^(/x(a;; — a)|g), let us 
consider the transformation formula 
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iaq, —iaq 



\ I o-2fc _„2 2fc+l -2n 

— a^g^ 



Q,-q ] = 391 



q',q'-+'], (5.29) 



which is true for any nonnegative integer values of A;. It is obtained by equating two expressions 
(5.22) and (5.25) for the dual discrete g-ultrasphcrical polynomials Dr^\fi(2k; —a)\q). Observe 
that (5.29) is still valid if one replaces numerator parameters q^'^^ and — a^g^'^+^ in both sides 
of it by c~^q~^'' and —ca^q'^'^'^^, c G C, respectively. Indeed, the left side of (5.29) represents 
a finite sum: 

In the case in question a = q~'^^ and (3 = — a^^^*^"*"^, so the g-shifted factorial {a,P;q)m in 
(5.30) is equal to 

{q-'\-aV'^';q)m 
m—l TO— 1 

= n [1 - - q'iq-''^ - a\^^'^^)] =J{[1- a'g'^+' - gV(2A;; -a^)], (5.31) 

3=0 j=0 

where, as before, fi{2k; —a^) = q~'^^ — a^q^^^^. The left side in (5.29) thus represents a 
polynomial Pn{x) in the fi{2k; —a?) of degree n. In a similar manner, one easily verifies that 
the right side of (5.29) also represents a polynomial p'n{x) in the same iJL{2k\—o?) of degree 
n. In other words, the transformation formula (5.29) states that the polynomials Pn{x) and 
Pn{x) arc equal to each other on the infinite set of distinct points Xk = fi{2k; —a^). Thus, 
they are identical. 

An immediate consequence of this statement is that (5.29) still holds if one replaces the 
numerator parameters and a^q'^'^'^^ in both sides of (5.29) by c~^q~'^^ and (x?q^^^^, 
respectively. The point is that 

m~l 
j=0 

m—l 

= ll[l-aV'^'-q^^i,{2k;-a% 

j=0 

where /ic(2fc; — a^) = c~^q~'^^ —ca^(^^'^^ . So, the replacements q"'^^ c'^q~'^^ and (]?(^^'^^ 
ca^q^^^^ change only the variable: fi{2k; —a^) /ic(2A;; — a^), whereas all other parameters 
in both sides of (5.29) are unaltered. Thus, our statement is proved. 

We are now in a position to establish other orthogonality relations for the polynomials 
Dli\lJ,{x;—a)\q), distinct from those, obtained in subsection 5.4. To achieve this, we use the 
fact that the polynomials Dn\n{x\ —a)\q) at the points x = x^^^ := 2k — \n{.y/aq/d)/\nq are 
equal to 

(y"(4 '-«)^) = 3^2 y /- ■ r q,-q ], (5-32) 
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where nix^^^; -a) = ^/aq [d ^ q "^^ - dq^^). From (5.27) and (5.29) (with q '^^ and -a^q^'^^^ 
replaced by d'~^q~'^^ and —da^q^^^^, respectively) it then follows that 

Hence, from the orthogonality relations (5.28) one obtains infinite number of orthogonality 
relations for the polynomials Dn\^^{x^,—a)\q), which are parameterized by the same d as in 
(5.28). They are of the form 

x£>W(M(x(f); -a)|g)^(«)(/x(x(f); -a)|g) = j^^^^ Srs, (5.33) 

I Q (^;Q )r 

where ti = ^ q^ /a and t2 = \/q/a. 

There exist yet another connection of the polynomials D^\ii{x)\q) with the polynomials 
(5.27). In order to obtain it we consider the relation 
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g-2fc-l^_„2^2fc+2^ 

iaq, —iaq 



— a^g^ 



q,-q ] = q 3<?>i 



q',q'^-'\. (5.34) 



This relation is true for nonnegative integer values of k. However, it can be proved (in the 
same way as in the case of formula (5.29)) that it holds also if we replace g^'^ and q~'^^ by 
cq^'' and c''^q~'^^, respectively. 

As in the previous case, wc put x = x^ := 2k — ln{y^aq/d)/ Inq, that is, fi[x^^'^ ; —a) = 
y/aq (yd^^ q^^^ — dq^^). Then by means of formula (5.34) (with q~'^^ and q^^ replaced by 
d~^q~'^^ and dq"^^, respectively), we derive that 



We may apply to the polynomials Unix; ^/ q/a, -sj q^ /a\q^) the orthogonality relations (5.28) 

and obtain an infinite number of orthogonality relations for the polynomials D^\ii{y)\q) . 
However, they coincide with the orthogonality relations (5.33). 

It is important to know whether an orthogonality measure for polynomials is extremal. 
The extremality of the measures in (5.33) for the polynomials Dii\fic{x', ~a.)\q) depends on 
the extremality of the orthogonality measures in (5.28) for the polynomials (5.27). If some 
of the measures in (5.28) are extremal, then the corresponding measures in (5.33) are also 
extremal. 

6. DUALITY OF BIG g-LAGUERRE AND g-MEIXNER POLYNOMIALS 

6.1. Operators related to big g-Laguerre polynomials 

Let H = Hahe the Hilbert space of functions, introduced in subsection 3.1, which is fixed 
by a real number a, such that < a < In this section we are interested in the operator 

Afn = rn+lfn+l + ^n/n-l + dnfn, (6-1) 



40 



where 

r-„ = (-a^g'^+^jVVll - - aq^)il - bq^), = -abq^-^+\l + q) + q''+\a + ah + b), 

a is the same as above and 6 is a fixed real number, such that 6 < 0. 

Since q < \ the operator A is bounded. Therefore, one can close this operator and we 
assume in what follows that A is a closed (and consequently defined on the whole space H.) 
operator. Since A is symmetric, its closure is a self-adjoint operator. In the same way as in 
subsection 3.1, one readily proves that A is a Hilbert-Schmidt operator. Therefore, A has a 
discrete spectrum. 

We wish to find eigenfunctions ^a(^) of the operator A, A^x{x) = ^^x{x)- We set 

oo 
n=0 

Acting by A upon both sides of this relation, one derives that 

oo oc 

anW[rn+l fn+1 + Tji fn-1 + C^n /n] = A ^ an{X)fn- 

n=0 n=0 

Comparing coefficients of a fixed /„, one obtains a three-term recurrence relation for the 
a„(A): 

(A) + (A) + dnan{X) = Aa„(A). 

Making the substitution 

a„(A) = (-a6)-"/29-"("+3)/i f ij^S^h}^\''\'^^X) 

and using the explicit expressions for the r„ and cZ„, we derive the relation 

(1 - aq^+'){l - 6g"+i)a;+i(A) - abq^+'{l - q^a'^^.iX) + d„<(A) = Aa;(A), 

where, as before, < a < q~^ and 6 < 0. It coincides with the recurrence relation for the big 
g-Laguerre polynomials, which are defined as 

Pn{\;a,b;q) := 3(/)2(g"", 0, A; aq,bq; q,q) 

= {Q-''/b;qV 201(9"", aq/X; aq; q,X/b) (6.2) 
(see formula (3.11.3) in [KSw]), that is, a^(A) = P„(A;a, 6;g). Therefore, 

1 /2 

a„(A) = (-a6)-/2Q-"("+3)/^ (^1^^) ^' ^^"^^ 

Thus, the eigenfunctions of the operator A have the form 

n=0 ^ y.1,1}n J 

= ^(_5)-/2„-3n/4^-n(n+3)/4 (^9ig)n ^)l/2p^(A; «, fe; g)x^ 
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Since the spectrum of the operator A is discrete, only a countable set of these functions 
belongs to the Hilbert space H. This discrete set of functions determines a spectrum of A. 

In order to be able to find a spectrum of the operator A, we consider the linear operator 
q-'^o -v^rhich is diagonal in the basis {/n} and is given by the formula 

(see formula (3.1)). We have to find how the operator q~'^° acts upon the eigenfunctions Cx{x) 
of the operator A (which belong to the Hilbert space Ti). In order to do this one can use the 
g-difference equation 

g-'^(l - g")A2 p„(A) = 5(A) PniqX) - [B{X) + Di\)]pn{\) + D{X) Pn{q-'X) (6.5) 
for the big g-Laguerre polynomials Pn{X) = PniX; a, b; q), where 

B{X) = abq{l - A), D{X) = (A - aq){X ~ bq). 

Multiply both sides of (6.5) by knfn{x) and sum up over n, where kn are the coefficients of 
Pn{X\a,b;q) in the expression (6.3) for the coefficients a„(A). Taking into account formula 
(6.4) and the form of the operator q~''° in the basis {/„}, one obtains the relation 

{aqf/\-'°XHx{x) = B{X)^,x{x) - [B{X) + D{X) - X'']^^) + D{X)^g-,^{x), (6.6) 

where B{X) and D{X) are the same as in (6.5). 

6.2. Spectrum of A and orthogonality of big g-Laguerre polynomials 

The aim of this subsection is to find a basis in the Hilbert space TC, which consists of 
eigenfunctions of the operator ^ in a normalized form, and to derive explicitly the unitary 
matrix U, connecting this basis with the canonical basis /„, n = 0, 1, 2, • • •, in 7^. This matrix 
directly leads to the orthogonality relation for the big q'-Laguerre polynomials. 

Exactly as in section 4, one can show that for some value of A (which belongs to the 
spectrum) the last term on the right side of (6.6) has to vanish. There are two such values of 
A: X = aq and A = bq, which are the roots of the equation -D(A) = 0. Let us show that both 
of these points are spectral points of the operator A. Due to (6.2) we have 

(_lj„)n n{n-l)/2 
Pniaq;a,b;q) = {q-yb;q)-'-^ ^ 



Pn{bq;a,b; q) 



{bq; q)n 
(-ag)'*g"("-^)/2 



(«?; q)n 

Hence, for the scalar product {^aq{x) , ^aq{x)) we have the expression 

{aq,bq;q)n „2/ , x 

^^{-abYq<n..V.^q,q^/niaq-,a,b-,q) 

= f; (-6/a)","("-i)/2 = iMa;b;q,b/a) = ^^/^ < oo. (6.7) 

{bq;q)n{q;q)n {b;q)oo 

Similarly, for {Cbq{x) , ^bq{x)) one has the expression 

E°° {aq,bq;q)n p2/L n {a/b;q)oo^ 
7 ,x„ n(n+3)/2 ( V^n \bq\ a, b] q) = r < OO. (6.8 
(-a6)"g"^"+'^)/^ (g; q)n {a; q)oo 
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Thus, the values X = aq and X = bq are the spectral points of the operator A. 

Let us find other spectral points of the operator A. Setting \ = aq in (6.6), we see that 

the operator q~'-^° transforms Caq{x) into a linear combination of the functions ^aq^i^) ^.nd 
We have to show that Caq'^i^) belongs to the Hilbert space H, that is, that 

/c c \ {aQ,bq;q)„ 2 , ^^ 

(W, W) =g (_,^,)n,„(„+3)/2(,.,)/n(a^ ; a, 6; ,) < oo. 

In order to achieve this we note that since (ag^; q)k = {aq; q)k{^ — aq^^^)/{l — aq), we have 



Pn{aq ■,a,h;q) 

k=o 



1 - ag^'+^ {q "■;q)k{aq]q)k 
^ l-aq {aq;q)k{bq;q)k {q;q)k 



1 yr-^ {q "■,q)k{aq;q)k q^ ( u \ 

<--. > -7 ^7 ^ = {l-aq) Pn{aq;a,b;q). 

wq)k{bq;q)k {q;q)k 

Therefore, the series for {^aq"^ i iaq^) is majorized (up to the finite constant (1 — aq)~^) by the 
corresponding series for {(,aq,iaq)- Thus, £,aq'^{^') cigcnfunction of A and the point aq^ 

belongs to the spectrum of the operator A. Setting A = aq^ in (6.6) and acting similarly, one 
obtains that .^^^3 (x) is an eigenfunction of A and the point aq^ belongs to the spectrum of 
A. Repeating this procedure, one sees that iaq-^{x), n = 1,2, • • •, are eigenfunctions of A and 
the set aq^, n = 1, 2, • • •, belongs to the spectrum of A. Likewise, one concludes that ^6gn(x), 
n = 1, 2, • • •, are eigenfunctions of A and the set hq^, n = 1, 2, • • •, belongs to the spectrum of 
A. So far we do not know whether the operator A has other spectral points or not. In order 
to solve this problem we shall proceed as in subsection 3.2. 

The functions ^oq"(a^) a^nd £,hq^{x), n = 1,2, • • •, are linearly independent elements of the 
Hilbert space Ji. Suppose that aq^ and bq^, n = 1, 2, • • •, constitute the whole spectrum of 
the operator A. Then the set of functions ^a(j"(a^) and ibq'^{x)-, ^ = 1) 2, • • •, is a basis in the 
space Ti.. Introducing the notations H„ := Cagn+i(a;) and := ibqn.+i{x), n = 0, 1, 2, • • •, and 
taking into account the relation B{\) + D{\) — = abq{l + q) — Xq{ab + a + b), we find from 
(6.6) that 

= a-3/26g-2"-3/2(i _ ag"+i)S„+i - a-^/^q-^''-^/^[b{l + q) - q^'+^ab + a + 6)]S„ 

+ a-3/26g-2"-V2(i _ qn^(^l _ aq^/b)En-l 
for A = aq^^^, that is, for ^aq"+''-{x) = ^n{x)- Similarly, 

= aV2rig-2"-3/2(i _ 6g-+i)s;+i - aV26-ig-2-3/2 
x[l + q- a-\''+\ab + a + b)]E'^ + aV2;,-ig-2n-i/2(^ _ ^n^^^ _ i,qya)E'^_^ 

for A = 6g"'+^, that is, for C6g"+i(^) = '^n(^)- 

As we see, the matrix of the operator in the basis E!„ = ^aqn+i(x), = ^i,gn+i{x), 
n = 0, 1,2,---, is not symmetric, although in the initial basis /„, n = 0, 1,2,---, it was 
symmetric. The reason is that the matrix M := {{bmn)m,n=o (^mn)m,n=o) "^^^^ entries 

bmn-=am{aq"-), b'^^:= am{bq'^), m, n = 0, 1, 2, • • • , 

where am{dq^), d = a,b, are coefficients (6.3) in the expansion S,dq"{x) = J2m '^m{dq"') fn{x) 
(see (6.4)), is not unitary. It maps the basis {/„} into the basis {^oqn+i,^6gn+i} in the Hilbert 
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space 7i = TCa- The nonunitarity of the matrix M is equivalent to the statement that the basis 
S„ := ^agn+i (x), := ^bqn+i (x), n = 0, 1, 2, • • •, is not normahzed. In order to normalize it, 
we have to multiply H„ by appropriate numbers c„ and by numbers c^. Let S„ = c„H„, 
= c[jH„, n = 0, 1, 2, • • •, be a normalized basis. Then the operator q~'^° is symmetric in 
this basis and has the form 

= c;^iC„a-3/26g-2«-3/2 (1 _ ag"+l)S„+i - a-3/2^-2n-3/2 [^(1 ^ ^) _ + „ + 

x[l + q- a-\^+\ab + a + b)]^^ + c';\c;ai/26-ig-2-i/2(i _ _ be/a)^^_^. 

The symmetricity of the matrix of the operator q~'^° in the basis S^} means that 

that is, 



Cn / (l-ag-) / (l-6g") 



Thus, 



and 



c„_i y "(1 -(?")(! -ag"/6)' \/ " (1 - - 6(/-/« 

1/2 



1 



where C and C are some constants. 
Therefore, in the expansions 

iaq»{x) = H„(x) = ^ c„ am{aq'^) fm{x) = ^ /m(a;), (6.11) 

m m 

4g„(a;) = S„(x) = XI amibq"") fm{x) = ^ fm{x), (6.12) 

m m 

the matrix M := ((6mn)^,n=o (^mn)m,n=o) with entries 

1 

bmn = Cn a^iaq^ = C { ^j^^ T^^^^T^^) ^ Q-"-^""^'^/' Pm{aq-^\a, b; q), (6.13) 

V {o-Q/b, q; q)n {q; q)m {-abj'^J 

1 

Knn=<i am{bql=C' ( j^T^ , ^ q-"-^"-^'^" Pm{bq-^\a, b; q), (6.14) 

\{{bq/a,q;q)n [q; q)m {-ab}'^ J 

is unitary, provided that the constants C and C are appropriately chosen. In order to calculate 
these constants, one can use the relations 



oo 



m=0 m=0 
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for n = 0. Then these sums are multiples of the sums in (6.7) and (6.8), so we find that 

{o/a;q)oi [a/b\q)^ 
The coefficients and in (6.11)-(6.14) are thus real and equal to 

\{aq/b,q]q)nib/a;q)oc,J 

and ^ 

^/ ^ / (^g;g)n(Qg;g)oog" \ ^ ^ 
" \ipq/a,q\q)n{a/h;q) 

The orthogonality of the matrix M = {bmn b'^n) means that 

bmnbmn' = Snn' , ^mn^mra' = ^nn' , '^bmnb'mn' = 0' (6.18) 

m m m 

^^{bmnbm'n + ^mn^m'n) = ^mm' ■ (6.19) 
n 

Substituting the expressions for bmn and into (6.19), one obtains the relation 

{bq;q)oo {aq;q)r 



E r..7'rr!-.^ Pm(ag"+^;a, 6; ,)P^, (ag-+^ a, 6; g) 



(6/a; g)oo ^ {aq/b; q)n{q; q)n 



. {aq;q)oo {bq;q)nq d ^i, n+i a n+i , N 

+ / \ 7n w ^Pm{bq ^ ■a,b;q)Pm'{bq ^ ■,a,b;q) 

{a/b-q)oo {bq/a;q)n{q;q)n 

= ( (-afe)"9"-^"+^)/Xm'. (6.20) 

[aq, bq; q)m 

This identity must give an orthogonality relation for the big g-Laguerre polynomials Pmiv) = 
Pm{y]0'^b;q). An only gap, which appears here, is the following. We have assumed that the 
points ag" and 6g", n = 0, 1, 2, • • •, exhaust the whole spectrum of the operator A. As in the 
case of the operator I2 in section 4, if the operator A had other spectral points Xfe, then on the 
left-hand side of (6.20) would appear other summands ^xp,Pm{xk] o,, b; q)Pm'{,Xk; a, b; g), which 
correspond to these additional points. Let us show that these additional summands do not 
appear. To this end we set m = m' = in the relation (6.20) with the additional summands. 
This results in the equality 

{bq;q)oo {aq;q)nq"- , {aq;q)oo {bq;q) 



V»g,gM'/ _^ l»gigjoo {vq,q)nq _^ y^ (6 21) 

{b/a;q)oo ^^{aq/b;q)n{q;q)n {a/b;q)oo '^^{bq/ a;q)n{q;q)n ^ ^'^ ' 

In terms of the 24'i basic hypergeometric series this identity can be written as 

{bq, g)oo (ag, 0; aq/b; g, g) + /"f ' 20i(^9, 0; bq/a; g, g) + V fix^, = 1- (6.22) 
(6/a;g)oo {a/b;q)oo ^ 

It represents a particular case of Sears' three-term transformation formula for 2^i{o-, 0; c; g, g) 
series (see [GR], formula (3.3.5)) if jixk = all values of k. Therefore, in (6.20) the sum 



45 



does really vanish and formula (6.20) gives an orthogonality relation for big g-Laguerre 
polynomials. 

By using the operators A and q~'''°, we thus derived the orthogonality relation for big 
g-Lagucrrc polynomials. 

The orthogonality relation (6.20) for big g-Laguerre polynomials enables one to formulate 
the following statement: The spectrum of the operator A coincides with the set of points aq^'^^ 
and bq^'^^, n = 0, 1, 2, • • •. The spectrum is simple and has one accumulation point at 0. 

6.3. Dual polynomials and functions 

The matrix AI = (bmn ^mn) with entries (6.13) and (6.14) is unitary and it connects two 
orthonormal bases in the Hilbert space TC. The relations (6.18) for its matrix elements is the 
orthogonality relation for the functions, which are dual to the big qi-Laguerre polynomials and 
are defined as 

/„(g-™;a,6|(?) :=P„(ag"+^a,6;Q), n = 0,l,2,---, (6.23) 

gn{q-'^;a,b\q):=Pm{bq''+^;a,b;q), n = 0,l,2,---. (6.24) 

Taking into account the expressions for the entries bmn and the first two relations in 
(6.18) can be written as 

oo oo 

^ a„(ag"+i)a„(ag"'+^) = c'^ d^n', Yl «m(6cz"+')a„.(5(/"'+i) = c'^' S^n'- 

m=0 m=0 

Substituting the explicit expressions for the coefficients am{Xn)-, we derive the following or- 
thogonality relations for the functions (6.23) and (6.24): 

E ^ ^7'^ ^2^^ g'"'^""^^^'gn(g-"';a,b|g)gn.(9-;a,fe|g) = c'-'Sr^r.', (6.26) 

m=Q ' 

E (,.g)!(!'at^)^ ^"^^"^"''^ f-^^'""-^ 3n' a, b\q) = 0, (6.27) 

where c„ and are given by the formulas (6.16) and (6.17). 
Comparing the expression 

M„(g--;a,6;g) := 2</.i(g-", g""; ag; g, -g"+ V?*) 

for the g-Meixner polynomials with the explicit form (6.2) of the big g-Laguerre polynomials 
Pm{x', a, b; g), we see that 

/„(g-">,6|g) = (g-™/6;g)-iM,(g-™;a,-fe/a;g). 

Since {q~"^/b; q)m = {bq; q)mi—b)~'^q~"^^'^~^^^^^, the orthogonality relation (6.25) leads to the 
orthogonality relation for the g-Meixner polynomials Mn{q~'^) = Mn{q~"^; a, —b/a; g): 

^0 {bq,q-q)m 
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= -77 ^ 7 ^ Q (>nn', [0.Z6) 

{oq; q)oo {aq\ q)n 

where, as before, < a < q"^ and 6 < 0. This orthogonahty relation coincides with the 
known formula for the g-Mcixner polynomials (see, for example, (3.13.2) in [KSw]). 

The functions (6.24) are also expressed in terms of g-Meixner polynomials. Indeed, we 
have 

5„(g-"';a,6|g) = 3</>2(g~'",0,6g"+^ aqM\ q,q) 

= (g-™/a;g)-S<^i(g-",g-"^; bq; qM^'^^/a) = {q-'^/a-q)-' Mn{q-^;b,-a/b; q), 

where 6 < 0, that is, one of the parameters in these g-Meixner polynomials is negative. 

Substituting this expression for gn{q~^;a,b\q) into (6.26), we obtain the orthogonality 
relation for g-Meixner polynomials Mn{q~"^) = Mn{q~'^;b, —a/b; q) with negative b: 

J. {bq-q)m{-a/br ^m{m-i)/2 M,{q-"^) M„,(g— ) 

{a/b;q)oo{hq/a,q-q)n ri;: /f. 0Q^ 

— q dnn'- (b.iyj 



(a«; 9)00 {bq; q)n 

Observe that this orthogonality relation is of the same form as for 6 > 0. As far as we 
know, this type of orthogonality relation for negative values of the parameter b has been first 
discussed in [AAK]. 

The relation (6.27) can be written as the equality 

°° /'_-|Nm m(m— 1)/2 

i \ M„(g— ;a,-6/a; q)M^,{q-^;b, -a/b; q) = 0, 

m=0 ^^"^ 

which holds for n,n' = 0, 1,2, The validity of this identity for arbitrary nonnegative 
integers n and n' can be verified directly by using Jackson's g-exponential function 

n=0 

and the fact that Eq{z) has zeroes at the points Zj = —q~-', j = 0,1,2, ■■■ . 

Notice that the appearance of the g-Meixner polynomials here as a dual family with 
respect to the big g-Laguerre polynomials is quite natural because the transformation q ^ q"^ 
interrelates these two sets of polynomials, that is, 

Mn{x; b, c; q'^) = {q'^-fb; q)n Pn{qx/b; 1/6, -c; q). 

Let us introduce the Hilbert space of functions F[q~"^) on the set m G {0, 1,2, • • •} 
with a scalar product given by the formula 

00 

{Fi,F2)b=J2 P(^) ^1 (^"") ^2(9-"^), (6.30) 

m=0 

where the weight function 

' {q\q)m{-abq^)^^ 
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is the same as in (6.25)-(6.27). Now we can formulate the following statement. 

Theorem 6.1. The functions (6.23) and (6.24) constitute an orthogonal basis in the 
Hilbert space ^. 

Proof. To show that the system of functions (6.23) and (6.24) constitutes a complete basis 
in the space ^ we take in the set of functions F^, k = 0,1,2,- ■ ■, such that Fj^{g~'^) = S^rn- 
It is clear that these functions constitute a basis in the space [^^. Let us show that each of 
these functions belongs to the closure V of the linear span V of the functions (6.23) and 
(6.24). This will prove the theorem. We consider the functions 



oo 



Fk{q-n = ^^knbmn + ^^'kn b'mn, A: = 0, 1, 2, • • • , 

n=0 n=0 

where bjn and 6^^ are the same as in (6.19). Then p~^(m)F}^{q~^) is an infinite linear 
combination of the functions (6.23) and (6.24). Moreover, Fk{q~'^) G V and, due to (6.19), 
Ffc, A; = 0, 1, 2, • • •, coincide, up to a constant, with the corresponding functions F^, introduced 
above. The theorem is proved. 

The weight function p(m) in (6.30) does not coincide with the orthogonality measure for 
g-Meixner polynomials. Multiplying this weight function by [{bq;q)m{—b)^'"'-q~^'^^^^'>/'^\~'^, 
we obtain the measure in (6.28). Let l^^^ be the Hilbert space of functions F{q~'^) on the set 
me {0, 1, 2, • • •} with the scalar product 

where the weight function coincides with the measure in (6.28). 

Taking into account the modification of the measure and the statement of Theorem 6.1, 
we conclude that the g-Meixner polynomials M„(g~'"; a, —b/a; q) and the functions 

{bq; q)m (-6)-™ ^^g— ; a, b\q) 

constitute an orthogonal basis in the space (^^^ . 

Proposition 6.1. The q-Meixner polynomials Mn{q~^;a,c;q), n = 0, 1,2,---, with the 
parameters a and c = —b/a do not constitute a complete basis in the Hilbert space i^^y that 
is, the q-Meixner polynomials are associated with the indeterminate moment problem and the 
measure in (6.28) is not an extremal measure for these polynomials. 

Proof. In order to prove this proposition we note that if the g-Meixner polynomials were 
associated with the determinate moment problem, then they would constitute a basis in the 
space of square integrable functions with respect to the measure from (6.28). However, this is 
not the case. By the definition of an extremal measure, if the measure in (6.28) were extremal, 
then again the set of the g-Meixner polynomials would be a basis in that space. Therefore, 
the measure is not extremal. Proposition is proved. 

Let now 1^2) be the Hilbert space of functions F(g~'") on the set m G {0, 1, 2, • • •}, with 
the scalar product 

(i^l,F2)(2) = f; q^rn^m F,{q--)F^). 
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The measure here coincides with the orthogonahty measure in (6.29) for g-Meixner polyno- 
mials Mn{q~"^;b,—a/b; q), b < 0. The following proposition is proved in the same way as 
Proposition 6.1. 

Proposition 6.2. The q-Meixner polynomials Mn{q~"^;b,—a/b;q), n = 0, 1,2, with 
b < do not constitute a complete basis in the Hilbert space l^g); that is, these q-Meixner 
polynomials are associated with the indeterminate moment problem and the measure in (6.29) 
is not an extremal measure for them. 

6.4. Generating function for big g-Laguerre polynomials 

The aim of this section is to derive a generating function for the big g-Laguerre polynomials 

G{x, t; a, b- q) := i^|^^g-(-i)/2p„(x; a, 6; q)t\ (6-31) 

n=0 ^'^ 

which will be used in the next subsection. Observe that this formula is a bit more general 
than each of the three instances of generating functions for big g-Laguerre polynomials, given 
in section 3.11 of [KSw]. 

Employing the explicit expression 

Pn{x;a,b;q) = {b~^q-'';q)~^2(l>i{<r'',aqx~^; aq; q,x/q) 

for big g-Laguerre polynomials, one obtains 

G{x,t;a,b; q) = y- :—{-bqt) > r I - 1 

t^o W<l)n ^ {aq,q;q)k ^bJ 

^ iaq,q;q)k{q;q)n-k 

{aq,q;q)k {q\q)m 
_ {aqx~^;q)k ^ ^^^k „-k{k-i)i2 \^ (ag*^+^;g)m ^ u^i-k+\m 

By the g-binomial theorem, the last sum equals to {—abq^;q)oo/{—bq^~''t;q)oo- Since 



{-bq'-''t;qU = q-^^^-^^'\-q/hqt;q)k{-hqt;q) 



OO) 



then 



Thus, 



-abq'^;q)^ {-abq^;q)^ q^'^^ 



{-bq^-H;q)oo {-bqt;q)^ {btf{-l/bt;q)k 



G{x, t; a, b; g) = — > ^ I - ) 

{-bqt;q)oo f~'{-l bt,q;q)k\bJ 



k=0 

,2.. 



L^^l^^cl^,(aqx-\ 0; -1/bt; q,x/b). (6.32) 
{-bqt; g)oo 
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This gives a desired generating function for big g-Laguerre polynomials. 
6.5. Biorthogonal systems of functions 

Note that the operator A from formula (6.1) can be written in the form 

where 

a={-ahqf'''{l-q), (3i = b{l+q), (32 = hq + aq{h + l), 
and J±, Q and q^° are the operators on H given as 

1/4 „/2 ^ 

= l_g V(l-g"+^)(l-«g"+^)/n+l, 

^-1/4 -(n-l)/2 ^ 

J- fn = V(l-9")(l-«9")/n-l, 

q'' fn = q'\aqf'^ fn, Qfn = q''fn. 
Prom the very beginning we could consider an operator 

Ai := aq-^^/^ [(1 - bQ)J+ + QJ-] - Piq'^'^° + P2Q, 

where a,(3i, and /?2 are the same as above. This operator is well defined, but it is not self- 
adjoint. Repeating the reasoning of section 3, we find that eigenfunctions of Ai are of the 
form 

n=n V W, q)n J 

= E a-^^'\-h)--l\-^^^^Pr.{X, a, b; (6.34) 

The last sum can be summed with the aid of formula (3.11.12) in [KSw]. We thus have 

M^) = {{-a/h'f^^x-q)^ ■ 2<^i(6gA-\ 0; hq- q,a-^''l\-hr^'\-^x\). 

Now we consider another operator 

A2 := aq^°/^ [J+Q + J-(l - hQ)] q^°/^ - Piq^^° + 

This operator is adjoint to the operator Ai : A2 = Ai. Repeating the reasoning of subsection 
6.1, we find that eigenfunctions of A2 have the form 

M^) = f;(-a6)-/\-"("+i)/2 fi^^^ifi^^)'^%„(A;a,6;,)/„(x) 

= f; a-3»/4(_b)-n/2g-n(n+l)/2 qUk; 'i)n p^^^. ^. ^^^n_ ^335) 
n=0 "^^^ 
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According to the formula (6.32), this function can be written as 

= fJ^^(-twt ^ ^M^Q/XA a^'\-hr'/\/x- q,X/b). 

[a ' [ 0) ' X, qjoo 

Let us denote by ^'^(a;), m = 0, ±1, ±2, • • •, the functions 

*m(a;) = c^'0aqm+i(x), m = 0, 1, 2, • • • , ■^_m{x) = c!^_-i^il}bqm{x), m = 1, 2, • • • , (6.36) 

and by (^m{x), m = 0, ±1, ±2, • • •, the functions 

^m{x) = Cm^aqm+i{x), m = 0, 1, 2, • • • , $_^(x) = c^_i(^6gm (x) , m = l,2,---, (6.37) 

where Cm and are given by formulas (6.9) and (6.10). 

Writing down the decompositions (6.34) and (6.35) for the functions ^m(a^) and $m(a^) 
(in terms of the orthonormal basis /„, n = 0, 1, 2, • • •, of the Hilbert space Ti) and taking into 
account the orthogonality relations (6.25)-(6.27) we find that 

{^m{.x), <^n{x)) = 6mn, m,n = 0, ±1, ±2, • • • . 

This means that we can formulate the following statement. 

Theorem 6.2. The set of functions ^'^(a;), m = 0, ibl,ib2, ■ • and the set of functions 
$m(x), m = 0, ±1, ±2, • • form biorthogonal sets of functions with respect to the scalar 
product in the Hilbert space 7i. 

7. ALTERNATIVE 9-CHARLIER POLYNOMIALS AND THEIR DUALS 

7.1. Pair of operators (i?i, J) 

Let Ti be the same separable complex Hilbert space as before. We have introduced into 
this space the orthonormal basis /„, n = 0, 1, 2, • • •, expressed in terms of monomials in x. We 
define on TL two operators. The first one, denoted as Q, acts on the basis elements as 

Qfn = g"/n. 

The second one, denoted as Bi , is given by the formula 

Blfn = Clnfn+l + ^n-l/n-l + bnfn, (7.1) 



with 



^ _ (a,^n+ui/2 V(l-g"+^)(TT^ 

^ ' (1 + ag2'^+i)^(l + ag2«)(l + ag2n+2) ' 



Xl+a92n)(l + ag2n+l) (1 + ag2n-l)(l + ag2n) 

where a is a fixed positive number. Clearly, Bi is a symmetric operator. 

Since a„ ^ and 6„ ^ when n — > oo, the operator Bi is bounded. Wc assume that it 
is defined on the whole Hilbert space 7i. For this reason, Bi is a self-adjoint operator. Let us 
show that Bi is a Hilbert-Schmidt operator. For the coefficients a„ and 6„ from (7.1), we have 
an+i/dn Q^^^ a^nd bn+i/bn — q when n oo. Since < ^ < 1, for the sum of all matrix 
elements of the operator Bi in the basis /„, n = 0, 1, 2, • • •, we have X^„(2a„ + 6„) < oo. This 
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means that 5 is a Hilbert-Schmidt operator. Thus, a spectrum of Bi is discrete and has a 
single accumulation point at 0. Moreover, a spectrum of Bi is simple, since Bi is representable 
by a Jacobi matrix with a„ 7^ (see [Ber], Chapter VII). 

To find eigenfunctions of the operator Bi, Bi^x = ^Cx^ we set 

n 

where /3n(A) are appropriate numerical coefficients. Acting by the operator Bi upon both 
sides of this relation, one derives that 

00 00 

XI /^'»(-^) («n/n+l + On-l/n-l + Kfu) = ^ XI /^"('^)/»' 
n=0 n=0 

where an and bn are the same as in (7.1). Collecting in this identity all factors, which multiply 
/„ with fixed n, one derives the recurrence relation for the coefficients PnW- 

/?„+i(A)o„ + /3„_i(A)a„_i +/3n(A)6n = A/?„(A). 

The substitution 

3 (X)- ( g)n (1 + V") \ -n(n+3)/4 g, (^-^ 

reduces this relation to the following one 

A _ n l+aq"^ _ 2n-l 1 ~ q"^ 



(l + ag2n)(l + aq2n+l)' -^'^ (1 + ag2n-l) (1 + ag2n) " 

This is the recurrence relation for the alternative g-Charlier polynomials 

Kn{\\a-q) :=2</>i(g"",-ag"; 0; q,q\) 
(see, formulas (3.22.1) and (3.22.2) in [KSw]). Therefore, /3;(A) = Kn{X;a;q) and 

PniX) = i ^T'^f^'^fn^ ) l-^'^^'^^'KniX; a; q). (7.2) 
V (g; g)n (1 + a)a" J 

For the eigenvectors we thus have the expression 

= X ( 7''Tnl7n ) a; «)/n. (7-3) 

^ V (9;?)n(H-a)a'* y 

Since the spectrum of the operator Bi is discrete, only for a discrete set of values of A these 
vectors belong to the Hilbert space H. 

Now we look for a spectrum of the operator Bi and for a set of polynomials, dual to 
alternative q'-Charlier polynomials. To this end we use the action of the operator 

J:=Q~^-aQ 
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upon the eigenvectors (,x, which belong to the Hilbert space H. In order to find how this 
operator acts upon these vectors, one can use the g-difference equation 



- ag")K„(A) = -aKniqX) + X'^KniX) - X'^l - X)Kn{q-'X) (7.4) 

for the alternative g-Charlier polynomials Kn{X) = Kn{X; a; q) (see formula (3.22.5) in [KSw]). 
Multiply both sides of (7.4) by fc„ and sum up over n, where kn are the coefficients of 
Kn{X;a;q) in the expression (7.2) for f3n{X). Taking into account the formula (7.3) and the 
fact that Jfn = {q~"' — aq'^)fn, one obtains the relation 

J^x = -a^gx + A-^Ca - A-^l - A) ^i;,, (7.5) 

which will be used in the next subsection. 



7.2. Spectrum of Bi and orthogonality of alternative g-Charlier polynomials 

The aim of this section is to find, by using the operators Bi and J, a basis in the Hilbert 
space Ti., which consists of eigenfunctions of the operator Bi in a normalized form, and to 
derive explicitly the unitary matrix U, connecting this basis with the basis fn, n = 0, 1, 2, • • •, 
in H. This matrix leads directly to the orthogonality relation for alternative g-Charlier poly- 
nomials. For this purpose we first find a spectrum of Bi. 

We proceed as in the previous cases. First we analyze a form of the spectrum of the 
operator Bi from the point of view of the spectral theory of Hilbert-Schmidt operators. If A 
is a spectral point of the operator Bi, then (as it is easy to see from (7.5)) a successive action 
by the operator J upon the function (eigcnfunction of Bi) ^\ leads to the eigenfunctions ^qm\, 
m = 0, ±1, ±2, • • •. However, since Bi is a Hilbert-Schmidt operator, not all of these points 
may belong to the spectrum of Bi, since g^™A — > oo when m — > +oo once A 7^ 0. This means 
that the coefficient 1 — A' of ^^-ia' in (7.5) must vanish for some eigenvalue A'. Clearly, it 
vanishes when A' = 1. Moreover, this is the only possibility for the coefficient of ^q-ix' in (7.5) 
to vanish, that is, the point A = 1 is a spectral point for the operator Bi. Let us show that 
the corresponding eigenfunction ^1 = ^^0 belongs to the Hilbert space 7i. 

By formula (II. 6) of Appendix II in [GR], one has Kn{l; a; q) = 24>i{q~"': —ciq""i 0; q, q) = 
(— a)"'g" . Therefore, 



i \ 9)n(l + T^2(-. x_Y^ (-a;g)n(l + Qg^")Q" 

^0 il+a){q-qU-qr^ir.+i)/2^n{^'^'1) A. il+q)iq;qU-ni3n-i)/2- 

In order to calculate this sum, we take the limit d, e ^ 00 in the equality 



(7.6) 



(1 + ag^")(~a; q)n{d; q)n{e; q)n raqy n{n-i)/2 ^ {-~aq\q)oo{-aq/de-q)^ 
^ (1 + a){-aq/d; q)n{-aq/e; q)n{q\ q)n ^deJ {-aq/d; g)oo(-ag/e; q)c 



n=0 

(see formula in Exercise 2.12, Chapter 2 of [GR]). Since 

lim (d;g),(e;(?),(ag/de)" = g"("-l)(ag)^ 

a,e— >oo 

we obtain from here that the sum in (7.6) is equal to {—aq', q)oo, that is, (^1,^1) < 00 and ^1 
belongs to the Hilbert space H. Thus, the point A = 1 does belong to the spectrum of the 
operator B^. 

Let us find other spectral points of the operator Bi (recall that a spectrum of Bi is 
discrete). Setting A = 1 in (7.5), we see that the operator J transforms into a linear 
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combination of the vectors and ^^o. Moreover, belongs to the Hilbert space H, since the 
series 

is majorized by the corresponding series (7.6) for ,^^0. Therefore, belongs to the Hilbert 
space H and the point q is an eigenvalue of the operator Bi. Similarly, setting A = g in 
(7.5), we find that is an eigenvector of Bi and the point belongs to the spectrum of Bi. 
Repeating this procedure, we find that all ^qn, n = 0, 1, 2, ■ ■ •, are eigenvectors of Bi and the 
set g", n = 0, 1, 2, • • •, belongs to the spectrum of Bi. So far, we do not know yet whether 
other spectral points exist or not. 

The functions ^qn, n = 0, 1,2, • • •, are linearly independent elements of the Hilbert space 
H (since they correspond to different eigenvalues of the self-adjoint operator Bi). Suppose 
that values g", n = 0, 1, 2, • • •, constitute a whole spectrum of Bi. Then the set of vectors 
^g", = 0, 1, 2, • • •, is a basis in the Hilbert space Ti. Introducing the notation := ^^fc, 
A; = 0, 1, 2, • • •, we find from (7.5) that 

JEk = -a Ek+i + g"^ Sfe - q-^{l - q^) 

As we see, the matrix of the operator J in the basis Sfc, A; = 0, 1,2, • • ■, is not symmetric, 
although in the initial basis |n), n = 0, 1,2,---, it was symmetric. The reason is that the 
matrix {amn) with entries a^n := Pmiq"'), m, n = 0, 1, 2, • • •, where Pmiq^) Sire the coefficients 
(7.2) in the expansion = Pm{q"')fn, is not unitary. This fact is equivalent to the 
statement that the basis E]„ = ^qn, n = 0, 1, 2, • • •, is not normalized. To normalize it, one has 
to multiply H„ by corresponding numbers c„. Let S„ = c„H„, n = 0, 1, 2, • • •, be a normalized 
basis. Then the matrix of the operator J is symmetric in this basis. Since J has in the basis 
{H„} the form 

JS„ = -c~l-^CnaEn+i + g"" S„ - c~i;^Cnq~"-{l - g") 

then its symmetricity means that c^^^c^a = c^^Cn-\-iq~"'~^{l — q"''^^), that is, Cn/cn-i = 
i/ag"/(l - g"). Therefore, 

c^ = c(aV('^+^)/V(9;5)n)^/^ 

where c is a constant. 
The expansions 

iqn{x) = En{x) = ^ Cn/3^(g")|m) = ^ amn\m) (7.7) 

m m 

connect two orthonormal bases in the Hilbert space Ti,. This means that the matrix {cLmn)-, 
m, n = 0, 1, 2, • • • , with entries 

amn = cMq ) = c ^ ^ ,^^.(^.1)/. j ^miq;a;q) (7.8) 

is unitary, provided that the constant c is appropriately chosen. In order to calculate this 
constant, we use the relation J2m=o l^mnP = 1 for n = 0. Then this sum is a multiple of the 
sum in (7.6) and, consequently, 

c= i-aq;q)^/'^. 
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The matrix (a^n) is real and orthogonal, that is, 

n m 

Substituting into the first sum over n in (7.9) the expressions for a^n, we obtain the identity 

°° „n^n(n + l)/2 ( r,n>>'-r-t\ n>>' ( r-f r-,\ 

E , a; q) a; ,) = "J/^"^^" ^-^-^^^/^ W , (7.10) 

^ (g;g)n (l + ag^'») 

which must yield the orthogonality relation for alternative g-Charlier polynomials. An only 
gap, which remains to be clarified, is the following. We have assumed that the points g", 

n = 0, 1, 2, • • •, exhaust the whole spectrum of Bi. Let us show that this is the case. 

Recall that the self-adjoint operator Ii is represented by a Jacobi matrix in the basis /„, 
n = 0, 1,2, • • •. According to the theory of operators of such type, eigenvectors ^\ of Bi are 
expanded into series in the basis /„, n = 0, 1, 2, • • •, with coefficients, which are polynomials in 
A. These polynomials are orthogonal with respect to some positive measure d^{X) (moreover, 
for self-adjoint operators this measure is unique). The set (a subset of R), on which these 
polynomials are orthogonal, coincides with the spectrum of the operator under consideration 
and the spectrum is simple. 

We have found that the spectrum of B\ contains the points g", = 0, 1, 2, • • •. If the 
operator Bi had other spectral points x, then on the left-hand side of (7.10) there would be 
other summands ^mixkl q) Kmi{xk', a; q), corresponding to these additional points. Let 
us show that these additional summands do not appear. We set m = m' = in the relation 
(7.10) with the additional summands. Since Ko{x;a;q) = 1, we have the equality 

-„n^n(n+l)/2 
n=0 y'i^'iJ'"' k 

According to the formula for the (/-exponential function Eq{a) (see formula (II.2) of Appendix 

II in [GR]), we have Xl^lio ""^"'^"''^^^^^(9; ^)n^ = i-<^QlQ)oo- Hence, J2kl^Xk = 0- This means 
that additional summands do not appear in (7.10) and it does represent the orthogonality 
relation for alternative g-Charlier polynomials. 

Due to the orthogonality relation for the alternative g-Charlier polynomials, we arrive at 

the following statement: 

Proposition 8.1. The spectrum of the operator Bi coincides with the set of points q^, 
n = 0, 1, 2, • • • . The spectrum is simple and has one accumulation point at 0. 

7.3. Dual alternative g-Charlier polynomials 

Now we consider the second identity in (7.9), which gives the orthogonality relation for 
the matrix elements cimn-, considered as functions of m. Up to multiplicative factors these 
functions coincide with 

Fn{x;a\q)=24>i{x,-a/x- 0; q,q''+^), (7.11) 
considered on the set x G {q~"^ | m = 0, 1, 2, • • •}. Consequently, 

iQ;i)n {-aq-^;qU{q;q)marr^q-^(rn-,i)/2 ^n[q ,a\q) 
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and the second identity in (7.9) gives the orthogonality relation for Fn{q '";a|g): 

The functions Fn{x;a,b\q) can be represented in another form. Indeed, taking in the 
relation (III. 8) of Appendix III in [GR] the limit c — ^ oo, one derives the relation 

20i(r™,-ar; 0; (Z, = (-arg-%0o(cZ-™,-or,r" -; 9,-«7«)- 

Therefore, we have 

Fniq-""; a\q) = {-arq'^'sMQ-'^, '^q^^, - ; Q, -97«)- (7-13) 

The basic hyper geometric function 3^0 in (7.13) is a polynomial of degree n in the variable 
fi{m) := q~"^ — aq"^, which represents a g-quadratic lattice; we denote it by 

d„(Mm);a;g) := 3<^o(g-™, -a g'", g""; -; g, -g7«) • (7-14) 
Then formula (7.12) yields the orthogonality relation 

^ (l + ag^"^)a"^ miSm-l)/2^^^ y ^^^^ y ^ = ^^iil^ J (7.15) 
^Q(-ag"»;g)oo(g;g)m ^ ^'^JnKI-K h Qn^n(n+l)/2 V J 

for the polynomials (7.14) when a > 0. We call the polynomials dn{lJi{m); a\ q) dual alternative 
q-Charlier polynomials. Thus, the following theorem holds. 

Theorem 7.1. The polynomials dn{l^{m); a; q) , given hy formula (7.14), are orthogonal 
on the set of points /x(m) := q~'^ — aq^, m = 0, 1,2, and the orthogonality relation is 
given by formula (7.15). 

Let be the Hilbert space of functions on the set m = 0, 1, 2, • • • with the scalar product 

= E ( ^^TiT\ q^^''^-'^"F,(m)F^), (7.16) 
(-ag'";g)oo(?;9)m 

where the weight function is taken from (7.15). The polynomials (7.14) are in one-to-one cor- 
respondence with the columns of the orthogonal matrix (d^^) and the orthogonality relation 
(7.15) is equivalent to the orthogonality of these columns. Due to (7.9) the columns of the ma- 
trix {cLrnn) form an orthonormal basis in the Hilbert space of sequences a = {a„ | n = 0, 1, 2, • • •} 
with the scalar product (a, a') = 'Yln^nO-'n- This scalar product is equivalent to the scalar 
product (7.16) for the polynomials dn{lJ-{m); a; q) . For this reason, the set of polynomials 
dn{lJ,{m);a;q), tt, = 0, 1, 2, • • •, form an orthogonal basis in the Hilbert space l^. This means 
that either the dual alternative q-Charlier polynomials dnip-im); a; q) correspond to determi- 
nate moment problem or the point measure in (7.15) is extremal if these polynomials corre- 
spond to indeterminate moment problem. This question will not be further pursued here. 

A recurrence relation for the polynomials dn{lJ.{m); a; q) is derived from (7.4). It has the 
form 

(g-"^ - aq^)dn{fi{m)) = -adn+Mm)) + g-"dn(M(m)) - g-"(l - g")dn-i(//(m)), (7.17) 
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where dn{lJ,{m)) = dn{fJ,{m);a;q). A g-difference equation for dn{lJ,{m); a; q) can be obtained 
from the three-term recurrence relation for alternative g-Charlier polynomials. 

Note that for the polynomials (i„(/i(m); a; with q < 1 we have the expression 

dn{^i{m)■,a■,q-')=3MQ~"',-aq'^,q-''■, 0,0; q,q). (7.18) 

However, the recurrence relation for these polynomials (which can be obtained from the 
relation (7.17)), does not satisfy the positivity condition A^Cn+i > 0, that is, they are not 
orthogonal polynomials for a > (as it is the case for alternative q'-Charlier polynomials). This 
positivity condition holds only if wc require a < 0. In this case, the polynomials (7.18) arc the 
continuous big g-Hermite polynomials Hn{x]a\q) (for an explicit form of these polynomials 
see, for example, [KSw], formula (3.18.1)), which are orthogonal on a certain continuous set. 

8. DUALITY OF AL-SALAM-CARLITZ I AND 
g-CHARLIER POLYNOMIALS 

8.1. Pair of operators {B2,Q~^) 

Let a be a real number such that a < 0. Let C be the separable complex Hilbert space 
with the orthonormal basis |n), n = 0,1,2,---. We define on C the operator B2, which is 
given by the formula 

B2\n) = a„|n + 1) + a„_i|n - 1) - bn\n), (8.1) 

with 

an = (-a)i/2g"/V(l - ^"+^), bn = {a + l)q^. 

Clearly, B2 is a bounded symmetric operator. Therefore, we assume that it is defined on the 
whole Hilbert space C For this reason, B2 is a self-adjoint operator. As in the previous cases, 
it is easy to show that B2 is a Hilbert-Schmidt operator. Thus, a spectrum of B2 is discrete 
and has a single accumulation point at 0. Moreover, a spectrum of B2 is simple, since B2 is 
representable by a Jacobi matrix with a„ 7^ 0. 

To find eigenfunctions $,x of the operator B2, B2^\ = \(,x, we set = J2n i^nWl'iT') , where 
/3„(A) are appropriate numerical coefficients. Acting by the operator B2 upon both sides of 
this relation, one derives that 

00 00 

/3n(A) [a„|n + 1) + an-i\n - 1) - bn\n)] = (3n{X)\n), 

n=0 n=0 

where a„ and 6„ are the same as in (8.1). Collecting in this identity all factors, which multiply 
|n) with fixed n, one derives the recurrence relation for the coefficients /3n{\)- 

Making the substitution 

/ -n{n-l)/2 

we reduce this relation to the following one 

/3;+i(A) + i-a)q^-\l - q")/?;_i(A) - (a + l)g"/?;(A) = A/3;(A). 
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This is the recurrence relation for the Al-Salam-CarUtz I polynomials 

UtHX;q) := (-a)-g"("-i)/22<^i(g-,A-i; 0; q,qX/a) (8.2) 
(see, formulas (3.24.1) in [KSw]). Therefore, = ut\\q) and 

/ -n(n-l)/2 

For the eigenvectors we thus have the expression 

Since the spectrum of the operator B2 is discrete, only for a discrete set of values of A these 
vectors belong to the Hilbert space £.. This discrete set of eigenvectors determines a spectrum 
of B2. 

Let us find a spectrum of the operator B2 and a set of polynomials, dual to Al-Salam- 
Carlitz I polynomials. For this purpose we use the operator Q~^, which is diagonal in the 
basis {|n)}, and is given as 

Q-i|n) = 

We have to find how the operator acts upon the eigenvectors ^x, which belong to the 
Hilbert space JC. To this end, one can use the g-difference equation for Al-Salam-Carlitz I 
polynomials, which can be written as 

q-''Ut\X;q)=aq-'X-^UtHqX-,q)-dxUt\X-q) + aX-\l-X)il-X/a)UtHQ-'>^-,q), (8-5) 

where dx = aA~^(l + - X - X/q). 

Multiply both sides of (8.5) by /c„ |n) and sum up over n, where /c„ are the coefficients of 

ut\X]q) in the expression (8.3) for /3n(A). Taking into account formula (8.4) and the fact 
that Q~^\n) = q~^\n), one obtains the relation 

Q-^ = aq-^X-%x - dxCx + aA-2(l - A/a)(l - A) ^i;,, (8.6) 

which is used in the next subsection. 

8.2. Spectrum of B2 and orthogonality of Al-Salam— Carlitz polynomials 

Let us analyze a form of the spectrum of B2. If A is a spectral point of the operator B2, 

then (as it is easy to sec from (8.6)) a successive action by the operator upon the vector 
(eigenvector of B2) ^x leads to the eigenvectors $,qmx, m = 0, ±1, ±2, • • •. However, since B2 is 
a Hilbert-Schmidt operator, not all of these points may belong to the spectrum of B2, since 
g~"^A — > 00 when m — +00 if A 7^ 0. This means that the coefficient of ^^-i^' in (8.6) must 
vanish for some eigenvalue A'. There are two such values of A: A = 1 and X = a. Let us show 
that both of these points are spectral points of i?2- Observe that [4,"'*(1; q) = (— a)"(3'"("^-'^)/^ 
and U^\a;q) = (— l)"-g"-("-~^)/2. Hence, for the scalar product (^1,^1) we have the expression 

°° -n(n-l)/2 °° -n(n-l)/2 
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Similarly, for (^aiCa) one has the expression 



(^"'^«) = E \n-n\ = (l/a;?)oo. (8.8) 

n=0 

Thus, the values A = 1 and A = a are spectral points of the operator B2. 

Let us find other spectral points of i?2- Setting A = 1 in (8.6), we see that the operator 
transforms ^^o into a linear combination of the vectors and ^1. We have to show that 
belongs to the Hilbert space C, that is, that 

00 n(n— 1/2 

It is made in the same way as in the case of the scalar product (V'aq^j'^aq^) in subsection 4.2. 
Therefore, belongs to the Hilbert space C and the point q is an eigenvalue of the operator 
-B2. Similarly, setting A = g in (8.6), we find that ^^2 is an eigenvector of B2 and the point 
belongs to the spectrum of B2. Repeating this procedure, we find that all ^^n, n = 0, 1, 2, • • •, 
are eigenvectors of B2 and the set gi", n = 0, 1, 2, • • •, belongs to the spectrum of B2. Likewise, 
one conclude that the elements ^aq", n = 0, 1, 2, • • •, are eigenvectors of B2 and the set aq^, 
n = 0, 1, 2, • • •, belongs to the spectrum of B2. So far, we do not know yet whether other 
spectral points exist or not. 

The vectors and ^aq", n = 0,1,2,- ■ ■, are linearly independent elements of the Hilbert 
space C Suppose that values q"^ and aq^, n = 0, 1,2,---, constitute a whole spectrum of 
B2. Then the set of vectors ^q-n. and ^ag", n = 0, 1, 2, • • •, is a basis in the Hilbert space C 
Introducing the notations := £,qk and S'^ := ^^^fe, /c = 0, 1, 2, • • •, we find from (8.6) that 

S„ = -dn^n + ag-'"(l - g")(l - g7«) 2n-l, 

where 

As we see, the matrix of the operator in the basis S„, S^, n = 0, 1, 2, • • •, is not symmetric, 
although in the initial basis |m), m = 0,1,2,---, it was symmetric. The reason is that the 
matrix M := ((a^„)^„^o (amn)m,n=o) = ((«mn) {a'^n)) with entries 

amn-= Pm{q^), a'^^ := Prn{aq^) , m,n = 0, 1,2, 

where (3rn{dq^), d = l,a, are the coefficients (8.3) in the expansion ^^g" = Pm{dq"')\'n) , 
is not unitary. This fact is equivalent to the statement that the basis H„ = ^qn, = ^aq"; 
n = 0,l,2,---, is not normalized. To normalize it, one has to multiply H„ by corresponding 
numbers c„ and by corresponding numbers c^. Let S„ = c„S„ and = c^S^, n = 
0, 1, 2, ■ ■ ■, be a normalized basis. Then the matrix of the operator is symmetric in this 
basis. Since has in the basis E'^} the form 

Q-^ En = c;;|ic„ag"^"-^ - dn H„ + c^iiC„ag"2"(l - g")(l - g"/a) E^-i, 
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then its symmetricity means that 

c-l,Cnaq-^--' = c-'cn+iaq-"il - - q^+'/a), 

c'-lic'^a-'q-'^-' = c'-'c'^^,a-\-'^-\l - aq^+'){l - q-+') 

that is, 



Cn-i V <-i \j {l-q^){l-aq^) 

Therefore, for the coefficients c„ and we have the expressions 



{Q/a;q)n'^{q;q)n'^ " (a?; 9)n''^(5; 9)^^ 

where c and c' are some constant. 
Thus, in the expansions 

iqn = S„ = ^ c„/3„(q")|m) = ^ a„„|m), (8.9) 

m m 

Lq^ = ^'n = Y^ c;/3^(ag")|m) = a'^nW)^ (8-10) 

m m 

the matrix M := {dmn <^'mn) with entries 

(„n „-m(m-l)/2 \ 

, ; ^ ^/ ^ ^ ,m ] U!;^i\q^-q\ (8.11) 

(^/a; 9)n(g; q)n) (q; q)m {-ar J 

(n -m{m-l)/2 \ ^^"^ 

7 w ^ w\ ^ ^m ] (8-12) 
[aq; q)n{q; q)n) {t, q)m {-ar J 

is unitary, provided that the constants c and c' are appropriately chosen. In order to calculate 
these constants, we use the relations X]^=o I'^mnP = 1 ^'^d Ylm=Q Wmnl"^ = 1 for n = 0. Then 
these sums are multiples of the sums in (8.7) and (8.8), so we find that 

c=(a;g)-V2, c' = (1/a; g)-V2. (8.13) 
The coefficients c„ and are thus real and equal to 

qn/2 ^ „n/2 



{a;q)lL^{q/a;q)l!'^{q]q)l!'^ " (1/a; 5)^^ (ag; 9)^^(5; g)^^ 
The orthogonality of the matrix M means that 

^ ^ a"innO"mn' ~ ^nn' 1 ^ ^ a"mnamn' ~ ^nn' j ^mn^^^/ = 0, (8.14) 

m m m 

y~^(QmnQm'n + ^mn^m'n) = '^mm'- (8.15) 
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Substituting into (8.15) the expressions for Oj^^j and a^„, we obtain 



1 OO 



\q;qW^''-^^/^Sr,n', (8.16) 



which must yield the orthogonality relation for Al-Salam-Carlitz I polynomials. A problem, 
which remains to be clarified, is the following. We have assumed that the points and ag", 
n = 0, 1, 2, • • •, exhaust the whole spectrum of B2- Let us show that this is the case. 

If the operator B2 had other spectral points Xk, then on the left-hand side of (8.16) there 

would be other summands fix,. [xk] q) {x^; q), corresponding to these additional points. 
Let us show that these additional summands do not appear. For this we set n = n' = in 
the relation (8.16) with the additional summands. This results in the equality 

°° 1 qin 



r 7 — ; r r h 7— r 7 r r h flxi. = 1- (8-17) 

(a; q)oo [q/a; q)m{q; q)m (1/a; q)oo [aq; q)m{q; q)m ^ 

In order to show that ^^x^^ = 0, take into account the relation 

{Aq/C,Bq/C;q)oo , , ^ rj n , {A,B-q)^ a ( a in n in 2 m \ 1 

——————— 2(/)i(^, B; C; q, q) + — — . „ , „ . 2MAq C, Bq C; q C; q,q) = l 

{q/C,ABq/C;q}oo {C/q,ABq/C;q)oo 

(see formula (2.10.13) in [GR]). Putting here A = 0, B = and C = q/a, we obtain relation 
(8.17) without the summand '^Zkl^Xk- Therefore, one has "^^.^J-Xk — 0- This means that 
additional summands do not appear in (8.16) and it does represent the orthogonality relation 
for the Al-Salam-Carlitz polynomials. Due to this orthogonality, we arrive at the following 
statement: 

Proposition 8.1. The spectrum of the operator B2 coincides with the set of points q^ 
and aq^, n = 0, 1, 2, • • •. This spectrum is simple and has one accumulation point at 0. 

8.3. Duals to Al-Salam— Carlitz I polynomials 

Now we consider the identities (8.14), which give the orthogonality relations for the matrix 
elements cimn and a^„, considered as functions of m. Up to multiplicative factors they coincide 
with 

Fn{x; a; q) = 2<t>i{x, g""; 0; q, g"+Va), n = 0, 1, 2, • • • , (8.18) 
F'^{x- a- q) = 2<^i(x, ag""; 0; q, q^+^),n = 0, 1, 2, • • • , (8.19) 
considered on the set of points q~'^, m = 0, 1, 2, • • •. Namely, we have 

(^n „— m(m— 1)/2 \ ' 

/ „n „-m(m-l)/2 \ 

&'mn = c' 7 ,\ ^ w w (-a)"^g"^("^-^)/2F;(g— ; a; g), 

V N; q)n{q; q)n) {q; q)m (-a)"^ / 
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where c and c' are given by (8.13). The relations (8.14) lead to the following orthogonality 
relations for the functions (8.18) and (8.19): 



(a; Q)J Y1 P(^)^n(.q a; q)Fn'{q a; q) = {q/a; q)n{q\ q)nq "^nn', (8.20) 
m=0 

oo 

(lA;^)oo p{m)F:,{q-'^;a-q)Fl^,{q-^';a;q) = iaq;qUq;q)nq-''Snn', (8.21) 

oo 

^ p(m)F„(g— ; a; g)F;, (g"™; a; g) = 0, (8.22) 

(_(j)"igm(m-l)/2 



m=0 



m=0 

where 

/9(m) 



(g;?)m 

Comparing the expression (8.18) for the functions Fn{q~"^;a;q) with the expression 

C„(g— ;a';g) := 2</'i(g-", g""*; 0; g, -g"+V«') 

for the g-Charlier polynomials, wc conclude that 

Fniq-""; a; q) = C„(g— ; -a; q). (8.23) 

Applying to the expression for the functions Fl^{q'~'^ ; a; q) the transformation formula 

20l(g-^fe; 0; g, z) = (fez/g)«20i(g-, g/z; 0; g,g/6) 

(see (III.6) from Appendix III in [GR]), we derive that 

FLiq"^; a: q) = {-ar^Cniq-"^; -1/a; q). (8.24) 

Substituting the expressions (8.23) and (8.24) into the relations (8.20) and (8.21), wc obtain 
the orthogonality relations for the g-Charlier polynomials Cn{q^"^; —a; q) and Cn{q^"^; — 1/a; g), 
where a < 0. For Cniq"^) = Cniq"^; a'; g), a' > 0, it has the form 

oo /m m{m—l)/2 

E ^T-^ Cn{q-nCn'{q-n = (-a';g)oog-"(-g/a';g)n(g;g)„<5„„,. 

m=0 ^^"^ 

It coincides with the orthogonality relation known from the literature (see, for example, Chap- 
ter 7 in [GR]). 

Thus, we have shown that duals of the family of Al-Salam-Carlitz I polynomials ujC^ (g~™; g) 
are two sets of g-Charlier polynomials, one taken with the parameter —a and the second one 
with the parameter —1/a. 

The relation (8.22) leads to the following equality for g-Charlier polynomials: 

<^ m(m-l)/2 

Yl VT^C'„(g-"^;-a;g)C,.(5"'";-V«;5)) = 0, 



m=0 



where a < 0. 
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The set of functions (8.18) and (8.19) form an orthonormal basis in the Hilbert space P 
of functions, defined on the set of points m = 0, 1, 2, • • •, with the scalar product 



(/i,/2) = p(m)/i(m)/2(m), 



m=0 

where p(m) is the same as in formulas (8.20)-(8.22). One can deduce from this fact that the 
q-Charlier polynomials Cn{cr"^',a!;q), a' > 0, correspond to indeterminate moment problem 
and the orthogonality measure for them, obtained above, is not extremal. 

9. DUALITY OF LITTLE g-LAGUERRE AND 
Al-SALAM-CARLITZ II POLYNOMIALS 

9.1. Pair of operators {B3,Q~^) 

Let H = Ha be a separable complex Hilbert space of functions, used in sections 3-7, with 
the polynomial basis /„, n = 0, 1, 2, • • •, in it. We fix a real number a such that < a < q~^. 
Let B3 be the operator on 7^ = Ha, acting upon the basis elements fn as 

Bsfn = fln/n+l + Qn-lfn-l + &n/n, (9-1) 

with 

a„ = -aV2g"+V2y(i _ ^n+i)(i _ 6„ = q^{l + a) - aq^''{l + q). 

Clearly, is a symmetric operator. 

Since a„ ^ and 6„ ^ when n ^ oo, the operator B^ is bounded. We assume that it 
is defined on the whole Hilbert space H and, therefore, it is a self-adjoint operator. Exactly 
as in the previous cases one can show that B^ is a Hilbert-Schmidt operator. This means 
that the spectrum of S3 is discrete and has a single accumulation point at 0. Moreover, a 
spectrum of i?3 is simple, since B^ is representable by a Jacobi matrix with a„ 7^ 0. 

To find eigenfunctions of the operator i?3, -63.^^ = ^Caj we set 

n 

where Pn{^) are appropriate numerical coefficients. Acting by the operator B^ upon both 
sides of this relation, one derives that 

00 00 

^ /3n(A) (On/n+l + ttn-lfn-l + bnfn) = PnWfn, 
n=0 n=0 

where a„ and 6„ are the same as in (9.1). Collecting in this identity all factors, which multiply 
/„ with fixed n, one derives the recurrence relation for the coefficients l3n{X): 

/3„+i(A)an + ^n-i(A)an-i + ^n(A)6n = A/3„(A). 

The substitution 

reduces this relation to the following one 

-g"(l - ag"+i)/3;+i(A) - aq^{l - q^Pn-iW + (?" - 0?'"+' + aq^ - aq^^)(3'^{\) = XP'^{X). 
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This is the recurrence relation for the Httle g-Laguerre (Wall) polynomials 

Pn{X;a\q) :=2</'i(g"'',0; aq; q;qX) = {a-^q-'';q)-^ 2</>o(9"'', -; q;\/a). (9.2) 
Thus, we have /?j[.(A) =Pn(A;a|g) and, consequently. 



This means that eigenfunctions of the operator are of the form 

iaq;q)k 

I 



The expression for the eigenfunctions can be summed up. To show this one needs to know a 
generating function 

F{x; t; a\q) := ^-^^ pn{x; a\q)e (9.5) 

n=0 

for the little g-Laguerre polynomials. To evaluate (9.5), we start with the second expression 
in (9.2) in terms of the basic hyper geometric series 20o- Substituting it into (9.5) and using 
the relation 



{q;q)k {q;q)k{q;q)n-k' 

one obtains that 

Interchanging the order of summations in (9.6) leads to the desired expression 

F{x;t;a\q) = Eq{-aqt)24)o{x~'^,Q; -; q\xt), (9.7) 

where Eq{z) = {—z;q)oo is the g-exponential function of Jackson. 

Similarly, if one substitutes into (9.5) the explicit form of the little g-Laguerre polynomials 
in terms of 2^1 from (9.2), this yields an expression 

Hx; t; a\q) = ^^7^ 2</'i(0, 0; q/t; q; qx). (9.8) 

Using in (9.4) the explicit form of the generating function (9.7) for the little g-Laguerre 
polynomials, we arrive at 

Ux)=Eq{-qa-^/^x)2M^~\ 0; -; g; a-^'^Xx). 

Another expression for i\{x) can be written by using formula (9.8). 

Since the spectrum of the operator i?3 is discrete, only for a discrete set of values of A the 
functions (9.4) belong to the Hilbert space H. This discrete set of eigenvectors determines a 
spectrum of B^. 

Now we look for a spectrum of S3 and for a set of polynomials, dual to little g-Laguerre 
polynomials. To this end we use the operator , where Q is given by the formula Qfn = 
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q^fn- In order to find how the operator Q ^ acts upon eigenfunctions of B^, one can use the 
g-difference equation 

g-"Ap„(A) = -apn{qX) + (2 + a - AK(A) - (1 - AK(g-U) (9.9) 

for the little (7-Laguerre polynomials Pn(A) = Pn(A;a|(7) (see formula (3.20.4) in [KSw]). 
Multiply both sides of (9.9) by d„ |n) and sum up over n, where d„ are the coefficients of 
Pn{\;a\q) in the expression (9.3) for (5n{\)- Taking into account formula (9.9) and the fact 
that Q~^fn = q~^fn, one obtains the relation 

Q-i 6 = -«A^' + \-\2 + a- A) - ^'Hl - A) ^^A, (9-10) 
which is used in the next section. 



9.2. The spectrum of and orthogonality of little q'-Laguerre polynomials 

Let us find, by using the operators i?3 and Q^^, a basis in the Hilbcrt space Ti., which 
consists of eigenfunctions of the operator B^ in a normalized form, and the unitary matrix A, 
connecting this basis with the initial basis /„, n = 0, 1, 2, ■ • •, in H. First we have to find a 
spectrum of B3. 

Let us first look at a form of the spectrum of B^. If A is a spectral point of the operator B3, 
then (as it is easy to see from (9.10)) a successive action by the operator upon the function 
(eigenfunction of S3) leads to the eigenfunctions ^qmx, m = 0, ibl,ib2, • • •. However, since 

is a Hilbert-Schmidt operator, not all of these points belong to the spectrum of B3, since 
q~"^\ 00 when m —>■ +00. This means that the coefficient A'~^ — 1 of Cq-^x' in (9.10) must 
vanish for some eigenvalue A'. Clearly, it vanishes when A' = 1. Moreover, this is the only 
possibility for the coefficient of ^^-i^' in (9.10) to vanish, that is, the point A = 1 is a spectral 
point for the operator B3. Let us show that the corresponding eigenfunction ^1 = ^^o belongs 
to the Hilbert space TC. 

One has the following equality 

(ag)"q"("-i)/2 

Pn(i;ak) = 2MQ ">0; a?; q^q) = — ? — ^ • 

[aq; Q)n 

Therefore, for the scalar product in H we have 

Thus, the point A = 1 does belong to the spectrum of B3. 

Let us find other spectral points of the operator B^. Setting A = 1 in (9.10), we see that 
the operator transforms ^^o into a linear combination of the vectors and ^^0. Moreover, 
belongs to the Hilbert space H, since the series 

is majorized by the corresponding series for ^^o. Therefore, belongs to the Hilbcrt space 
H and the point q is an eigenvalue of the operator B^. Similarly, setting \ = q in (9.10), one 
finds likewise that is an eigenvector of B^ and the point q^ belongs to the spectrum of ^3. 
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Repeating this procedure, we find that all ^^n, n = 0, 1, 2, • • •, are eigenvectors of and the 
set g", n = 0, 1, 2, • • •, belongs to the spectrum of B^. So far, we do not know yet whether 
other spectral points exist or not. 

The functions ^^n, n = 0,1,2,---, are linearly independent elements of the Hilbert space 
7i. Suppose that values g", n = 0, 1, 2, • • •, constitute a whole spectrum of B^. Then the set 
of functions ^qn, n = 0, 1,2, • • •, is a basis in the Hilbert space H. Introducing the notation 
"fc ^qk, k = 0,1,2,- ■ ■, we find from (9.10) that 

Q-i Hfc = -aq-" Ek+i + g-'=(2 + a - g^) - g-'=(l - g'^) E^-i. (9.12) 

As we see, the matrix of the operator Q^^ in the basis Ej^, k = 0,1,2, is not symmetric, 
although in the initial basis |n), n = 0, 1,2, • • •, it was symmetric. The reason is that the 
matrix (amn) with entries Qmn '■= Pm{(t')) m,n = 0,1,2, - ■■ , where /?m(9") are the coefficients 
(9.3) in the expansion = Pm{Q^)fm is not unitary. This fact is equivalent to the 
statement that the basis H„ = ^qn, n = 0, 1, 2, • • •, is not normalized. To normalize it, one has 
to multiply H„ by corresponding numbers c„. Let S„ = c^S^, n = 0, 1, 2, • • •, be a normalized 
basis. Then the matrix of the operator is symmetric in this basis. It follows from (9.12) 
that has in the basis the form 

<5~^ S„ = -c~|iCng~"a + g-"(2 + a - g") S„ - c-^Cnq'^'il - g") 

The symmetricity of Q^^ in the basis means that c~^iCnaq~"' = c~^c„+ig~"~"^(l — g"^"^), 
that is, Cn/cn-i = y/aq/{l — g"). Therefore, 

Cn = ciiaqr/{q;q)n)'/^ 

where c is a constant. 
The expansions 

^qn(x) = ^fi(^X^ = ^ ^ CjiPuii^q )fm = ^ ^ O-rnnfm 
m m 

connect two orthonormal bases in the Hilbert space H. This means that the matrix (amn), 
m,n = 0,1,2, ■■■ , with entries 

is unitary, provided that the constant c is appropriately chosen. In order to calculate this 
constant, we use the relation Y1'^=q I^OnP = Yl'^=Q^nl^oi^"^) ~ ^- Since ^^{q^) = 1 and 

we have 

c= (ag;g)^^. 

The matrix A := (amn) is real and orthogonal. Thus, if S„, = 0, 1, 2, • • •, is a complete 
basis in H, then AA~^ = A~^A = E, that is, 

^ ^ amnO'm' n ~ ^mm' ^ ^ ^ O'mnO'ran' ~ ^nn' • ('^•-'-^) 
n m 
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Substituting into the first sum over n the expressions for CLmn, we obtain the identity 

f2j^Pm{Q'';a\q)pm'{q-;aq) = ^^^l^^f^ 5^^. , (9.15) 

[q; q)n [aq; qjooim; q)m 

which must yield the orthogonahty relation for little g-Laguerre polynomials. However, we 
have assumed that the points g", n = 0, 1, 2, • • •, exhaust the whole spectrum of B-,,. Let us 
show that this is the case. The reasoning here is the same as in the previous sections. Namely, 
we have found that the spectrum of S3 contains the points g", n = 0, 1, 2, • • •. If the operator 
had other spectral points x^, then on the left-hand side of (9.15) there would be other 
summands ^JixkPm{xk',o\q)p,m'{xk',o\q) with positive ^a;^., corresponding to these additional 
points. Let us show that these additional summands do not appear. We set m = m' = in 
the relation (9.15) with the additional summands. Since po{x]a\q) = 1, we have the equality 

n=0 k 

This formula is true if jjLx^. = 0. This means that additional summands do not appear in 
(9.15) and thus (9.15) does represent the orthogonality relation for little g-Laguerre polyno- 
mials. Consequently, the following proposition is true: 

Proposition 9.1. The spectrum of the operator B3 coincides with the set of points q^ , 
n = 0, 1, 2, • • • . This spectrum is simple and the functions ^q-n., n = 0, 1, 2, • • form a complete 
set of eigenfunctions of B3. The matrix (a^n) with entries (9.13) relates the initial basis {fn} 
with the normalized basis 

9.3. Al-Salam— Carlitz II polynomials as duals to little g-Laguerre polynomials 

Now we consider the second relation in (9.14). Taking into account the explicit expression 
for dmn, one obtains the orthogonality relation for the functions 

F„(g— ; a\q) := (fl-^g— ; q)'' 9""; - ; Q, 97«)- (9-16) 

This relation has the form 

f2^aq)-"^^^^F^iq-^; a\q)Fr,,iq-^; a\q) = {aq)-^S^^5nn'. (9.17) 

Comparing (9.16) with the Al-Salam-Carlitz II polynomials 

V;(«)(x;g) = (-a)"g-"("-i)/22<^o(g-", x; - ; g,g7a). 

we see that they are related to the functions (9.16), and (9.17) therefore leads to the orthog- 
onality relation for the Al-Salam-Carlitz II polynomials 



2 

E g"" T^(a) (-m. „^.T^(a) r-m. ^ _ ""(g; 9)« . 
7 ^ — 7 Vl 'gJV W 'gJ — 7 ^ ZJ'^nn') 
^^Q{(l\(l)m{aq\q)m " (ag; g)oog" 



known from the literature. Thus, Al-Salam-Carlitz II polynomials are duals to the little 
g-Laguerre polynomials. 
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Appendix 

In this appendix we prove the summation formula 



E 

n=0 



{abq, bq; q)n 1 - Q6g^"+^ 
{aq, q; q)n 1 - ahq 



{abq^;q)c 



(Al) 



First of all, observe that when 6 = this relation reduces to 



°° o"g"' 1 



n=0 



{aq, q; q)n {aq; q) 



which is a well-known limiting form of Jacobi's triple product identity (see [4], formula (1.6.3)). 
One can employ an easily verified relation 



(gg, -gg; q)n ^ 1 - a q 
(g, -g; g)„ 1 - g^ 



2„2n 



{A.2) 



in order to express the infinite sum in (A.l) in terms of a very-well-poised 4^5 basic hyperge- 
ometric series. This results in 



{abq, bq; q)n 1 - g6g2"+i „ ^2 
/ , r ^ ; zn — 1 = 4.(P5 



n=0 



{aq, g; q)n 1 - abq 



abq, bq, qy/abq, —qy/abq 
aq, \/abq, —\/abq, 0, 



g,gg 



(^.3) 



The next step is to utilize a limiting case of Jackson's sum of a terminating very- well-poised 
balanced ^4>^ series, 



6^5 



a, q^/a, —q\/a, b, c, d 
y/a, —\/a, aq/b, aq/c, aq/d 



aq 
bed 



{aq, aq/bc, aq/bd, aq/cd; g)oo 
{aq/b, aq/c, aq/d, aq/bcd; g)oo ' 



(^.4) 



which represents a g-analogue of Dougall's formula for a very- well-poised 2-balanced y-Fg series. 
When the parameters c and d tend to infinity, from (A. 4) it follows that 



a, qy/a, ~qy/a, b 
y/a, -y/a, aq/b, 0, 



{aq;q)c 



aq 

b ) {aq/b;q)c 



{A.5) 



To verify this, one needs only to use the limit relation 

With the substitutions a — > g6g and 6 — > 6g in (A.5), one recovers the desired identity (A.l). 
Similarly, when d — ^ 00 we derive from (A.4) the identities 



(1 - abq-" ^-){aq.. abq/c. abq: q)n n(n-i)/2 _ {abq^ ■ <'/"■ 'l) c 



„ (1-gi 



n=0 



(1 - abq){bq, cq, q; q)n{-a/cY 



{bq, cq; g)c 



(1 - abq^"-+^){abq, bq, cq; q)n n(n-i)/2 _ {abq^,a/c;> 



n=0 



(1 - abq) {aq, abq/ c, q; q)n{-c/ g)" 



{aq, abq/c; q) 



00 
00 



(^.6) 
(^.7) 



They have been employed in section 7. 
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We conclude this appendix with the following remark. There is another proof of the 
identity (4.33), based on vital use of the same summation formula (A.4). Actually, a relation 
may be derived, which is somewhat more general than (4.33). Indeed, consider the function 



%(a;,) := ^(-1) V^^^^/^ ' T 7'7 /(n;a) (A8) 

^j^Q ^q Wi 9)n 

for arbitrary nonnegative integers k, where the g-quadratic lattice /x(n; a) is defined as above: 

/x(n; a) := g-'* + ag"+\ (A9) 

We argue that all ryfe(a; <;) = 0, = 0, 1, 2, • • •. To verify that, begin with the case when A; = 
and employ relation (A.2) to show that 



770 (a; g) = 303 



QVaq, -q^jaq, aq 
aq, -V^, 



9,1 • 



The summation formula (A.4) in the limit as d ^ oo takes the form 

aq\ {aq,aq/bc;q)c 



595 



a, q^/a, -q^/a, b, c 
^/a, -\/a, aq/b, aq/c, 



be J {aq/b, aq/c; q)c 
In the particular case when be = aq this sum reduces to 



(^.10) 



393 



a, q^/a, —q\fa 
0, v^, —^Ja 



{aq,l;q)oo 

QA] = -77 ^ =0, 

(6, c; q)oo 



since {z;q)oo = for z = 1. Consequently, the function r)o{a;q) does vanish. 

For k = 1,2,3, one can proceed inductively. Employ the relation q^{n + l;a) = 
li{n\ q^a) to show that 

m+i{a;q) = (1 + aq)rik{a;q) - q~''~^{l - aq'^){l - aq^)r}k{aq^ \ q) ■ 

So, one obtains that indeed all rjk{a;q), k = 0, 1,2, vanish. The identity (4.33) is now 
an easy consequence of this statement if one takes into account that a product of the two 
polynomials D„(/x(m); a, 6, c|g) and Dn'{^x{m);b,a,abq/c\q) in (4.33) is some polynomial in 
/i(m) of degree n + n' . This completes the proof of (4.33), which is independent of the one, 
given in section 4. 
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